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This paper explains the multi-orbital band structures and itinerant magnetism of the iron-pnictide and 
chalcogenide superconductors. We first describe the generic band structure of a single, isolated FeAs layer. 
Use of its Abelian glide-mirror group allows us to reduce the primitive cell to one FeAs unit. For the 
lines and points of high symmetry in the corresponding large, square Brillouin zone, we specify how the 
one-electron Hamiltonian factorizes. From density-functional theory, and for the observed structure of 
LaOFeAs, we generate the set of eight Fe d and As p localized Wannier functions and their tight-binding 
(TB) Hamiltonian, h (k). For comparison, we generate the set of five Fe d Wannier orbitals. The topology 
of the bands, i.e. allowed and avoided crossings, specifically the origin of the d 6 pseudogap, is discussed, 
and the role of the As p orbitals and the elongation of the FeAs4 tetrahedron emphasized. We then couple 
the layers, mainly via interlayer hopping between As p z orbitals, and give the formalism for simple tetrag- 
onal and body-centered tetragonal (bet) stackings. This allows us to explain the material-specific 3D band 
structures, in particular the complicated ones of bet BaFe2 As2 and CaFe2As2 whose interlayer hoppings 
are large. Due to the high symmetry, several level inversions take place as functions of k z or pressure, and 
linear band dispersions (Dirac cones) are found at many places. The underlying symmetry elements are, 
however, easily broken by phonons or impurities, for instance, so that the Dirac points are not protected. 
Nor are they pinned to the Fermi level because the Fermi surface has several sheets. From the paramagnetic 
TB Hamiltonian, we form the band structures for spin spirals with wavevector q by coupling h (k) and 
h (k + q) . The band structure for stripe order is studied in detail as a function of the exchange potential, 
A, or moment, m, using Stoner theory. Gapping of the Fermi surface (FS) for small A requires matching of 
FS dimensions (nesting) and d-orbital characters. The interplay between pd hybridization and magnetism is 
discussed using simple 4x4 Hamiltonians. The origin of the propeller-shaped Fermi surface is explained 
in detail. Finally, we express the magnetic energy as the sum over band-structure energies and this enables 
us to understand to what extent the magnetic energies might be described by a Heisenberg Hamiltonian, 
and to address the much discussed interplay between the magnetic moment and the elongation of the FeAs4 
tetrahedron. 

Copyright line will be provided by the publisher 

['jMax-Planck-Institute for Solid State Research, Heisenbergstrasse 1, D70569 Stuttgart, Germany 



Foreword 

This paper is dedicated to Manuel Cardona, on the occasion of his 75th birthday. Manuel Cardona is world- 
famous not only for being an outstanding and creative experimentalist, but also for having deep theoretical 
insights, in particular into the band structures of solids. His book [1] with Peter Yu on Fundamentals of 
Semiconductors, now a classic of solid-state physics, devotes a large part to the explanation of theory in 
simple, accessible language. For our contribution to the present Festschrift, we have chosen a topic off the 
main line of this special issue, but we hope that Manuel and others will appreciate our effort to follow his 
example and make a complicated problem more transparent. 

We started working on superconductivity in the iron pnictides, and the later discovered chalchogenides, 
not only because we felt that this phenomenon would create the same kind of excitement as high-temperature 
superconductivity in cuprates had done two decades earlier, an excitement that the senior author had ex- 
perienced under the leadership of Manuel Cardona, but also because we -like many other researchers- 
had hoped that this would provide a shortcut to understanding high-temperature superconductivity. In fact, 
these new iron-based superconductors and the cuprates (as well as heavy-Fermion systems) share the most 
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important property of having a superconducting phase close to -or even coexisting with- an antiferromag- 
netic phase. In the iron-based superconductors, however, the latter is metallic rather than insulating, and 
these superconductors may therefore lack the complications of strong-correlation Mott-Hubbard physics. 
In both iron-based superconductors and cuprates, phonons seem to play a minor role for the coupling, 
although the debate is still alive @. After two and a half years with enormous theoretical and experi- 
mental effort, published in over 2000 papers, (31 S3 [5] it has become clear that, although superconductivity 
in the iron-based superconductors could turn out to be simpler than in the cuprates, the former are in 
many respects more complicated than the latter -their band structures, for example, are considerably more 
intricate- and our present understanding is far from complete. 

Trying to swim upwards through this cascade of papers has blinded at least the senior author and made 
him focus on finally publishing our own results. Encouraged by Manuel Cardona, to whom this could 
never have happened, we now take the opportunity to publish a pedagogical, self-contained account of the 
band structures and itinerant magnetism in the iron-based superconductors, stressing the role of symmetry 
and, as far as possible, reducing problems to that of diagonalizing a 2 x 2 matrix. 



1 Introduction 

The first report of superconductivity in an iron pnictide, specifically in F-doped LaOFeP below 5 K in 
2006 (6] 13, was hardly noticed and only two years later, when F-doped LaOFeAs was reported to su- 
perconduct below 28 K, the potential of iron pnictides as high-temperature superconducing materials was 
realized. [8| Following this discovery, more than 50 new iron superconductors with the same basic struc- 
ture were discovered (9] with T c reaching up to 56 K. iflOl This structure is shown in Fig. Q]for the case 
of LaOFeAs. The common motive is a planar FeAs layer in which the Fe atoms form a square lattice, 
tetrahedrally coordinated with As atoms placed alternatingly above and below the hollow centers of the 
squares. Instead of As, the ligand could be another pnictogen (P) or a chalcogen (X=Se or Te), but for 
simplicity, in this paper we shall refer to it as As. These superconductors are divided in four main families 
depending on their 3D crystal structure ifTTl : The iron chalcogenides are simple tetragonal (st) with the 
FeX layers stacked on top of each other (11 family). The iron pnictides have the FeAs layers separated by 
alkali metals (111 family), or by rare-earth oxygen/fluoride blocking layers (1111 family as in Fig.[T]i, in st 
stacking, or by alkali-earth metals (122 family) in body-centered tetragonal (bet) stacking. 

Iron-based superconductors share some general physical properties, although the details are often spe- 
cific to families, or even to compounds. With the exception of LiFeAs, the undoped compounds are 
spin-density wave (SDW) metals at low temperature with the Fe spins ordered anti-ferromagnetically be- 
tween nearest neighbors in the one direction and ferromagnetically in the other, thus forming stripe or 
double-stripe (FeTe) patterns. The values of the measured magnetic moments range from 0.4 ^ts/Fe in 
LaOFeAs 1121 . to ~ 1 fj, B in BaFe2As2 compounds, to over2 {1b in doped tellurides. J3]|4l|5| At a temper- 
ature above or at the Neel temperature, which is of order 100 K, there is a tetragonal-to-orthorhombic phase 
transition in which the in-plane lattice constant contracts by 0.5-1.0% in the direction of ferromagnetic or- 
der. Superconductivity sets in when the magnetic order is suppressed by pressure, electron or hole doping, 
or even isovalent doping on the As site, and at a much lower temperature. Both superconductivity and mag- 
netism are found to depend crucially on the details of the crystal structure; for example is it often observed 
that the highest T c s occur in those compounds where the FeAs4 tetrahedra are regular. BX3J Critical temper- 
atures range from a few K in iron-phosphides to 56 K in SmOFeAs. The variations in the phonon spectra 
are, however, small and seem uncorrected with T c . This, together with the proximity of magnetism and 
superconductivity in the phase diagram, was a first indication that the superconductivity is unconventional. 
A stronger indication seems to come from the symmetry of the superconducting gap, which is currently a 
strongly debated issue. lfl4l Depending on the sample, and on the experimental technique, multiple gaps 
with s symmetry and various degrees of anisotropy -but also of nodes- have been reported. J3]|4j|5l It now 
seems as if the gap symmetry is not universal, but material specific in these compounds. 
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Fig. 1 The layered structure of simple tetragonal LaOFeAs. The 3D primitive cell contains one Fe2As2 and one 
La2 02 layer, each containing three sheets: a square planar Fe (red) or O (blue) sheet sandwiched between two planar 
As (green) or La (yellow) sheets. c=874 pm. The coordination of Fe with As, or O with La, is tetrahedral. x and y 
are the vectors between the Fe-Fe or O-O nearest neighbors (separated by a=285 pm) and X and Y are those between 
As-As or La-La nearest neighbors in the same sheet. The directions of those vectors we shall denote x, y, X,and Y. 



Current understanding of the basic electronic structure has been reached mainly by angle-resolved pho- 
toemission (ARPES) |H5] [J6] [121 [J8] [J9] EH1 EI] E2] E3 |24l , quantum oscillation, and de-Haas-van- Alphen 
(dHvA) experiments 12511261 127 28 29][30][3T) in combination with density-functional (DFT) calculations. 
02 [33l [36l |37l EH1 E9l SOI Sll SZl S21 AU parent compounds have the electronic configuration Fe 

d 6 and are metallic. In all known cases, the Fermi surface (FS) in the paramagnetic tetragonal phase has 
two concentric hole pockets with dominant d xz /d yz character and two equivalent electron pockets with 
respectively d xz / d xy and d yz / d xy character. A third hole pocket may also be present, but its character, 
d xy or d 3z 2_ 1 , as well as the sizes and shapes of all sheets, vary among different families of compounds, 
and, within the same family, with chemical composition and pressure. In all stoichiometric compounds, 
the volumes of the hole sheets compensate those of the electron sheets. The magnetically stripe-ordered 
phase remains metallic, but the FS becomes much smaller and takes the shape of a propeller |fT8l plus, 
possibly, tiny pockets. J25) 

Given the strong tendency to magnetism, and the low value of the calculated electron-phonon interac- 
tion, l44l l45l 1461 spin fluctuations are the strongest candidate for mediating the superconductivity. Al- 
ternative scenarios have been proposed, in which superconductivity is due to magnetic interactions in the 
strong-coupling limit, polarons, or orbital fluctuations. l47l Models for spin fluctuations are based on the 
weak-coupling, itinerant limit, with superconductivity related to the presence of strong nesting between 
hole and electron sheets of the paramagnetic Fermi surface, which is also held responsible for the insta- 
bility towards magnetism. H2l l43l l48l This possibility has been investigated using more ore less sound 
models of the band structure, combined with different many-body methods (RPA, FLEX, frG, model ME 
calculations) which do seem to agree on a picture with competing instabilities towards magnetism and 
superconductivity. S3 |50l [3TI [32 [53 E^ [33 [5^ [53 [58l [5^ [60) The superconducting phase should 
be characterized by multiple gaps, with s and d symmetries almost degenerate. Modifying the shape and 
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orbital characters of the different sheets of the Fermi surface by doping, pressure, or chemistry can in- 
fluence the leading instability and affect the structure of the gap. As a result, a reasonable, qualitative 
picture of the materials trend, such as the dependence of T c and gap symmetry on the tetrahedral angle, 
has evolved. Il52l l57l Most experimental evidence seems to support this picture, but several points remain 
controversial. A badly understood issue is how to include 3D effects, which is particularly serious for the 
bet 122 compounds. 

Another problem concerns the magnetism: While it is true that spin-polarized DFT (SDFT) calculations 
reproduce the correct atomic coordinates and stripe-order of the moment, the magnitude of the moment is, 
except in doped FeTe, at least two times larger than what is measured by neutron scattering, or inferred 
from the gaps measured by ARPES, [30, 61|,dHvA, and optics, ll62l albeit much smaller than the satura- 
tion moment of 4 /is/Fe. Suppressing the too large moments in the calculations will, however, ruin the 
good agreement for the structure and the phonon spectra. l45l l63l [64 65 , 66j This over-estimation of the 
moment is opposite to what was found 25 years ago for the superconducting cuprates where the SDFT gave 
no moment, but is typical for itinerant magnets close to a magnetic quantum critical point (QCP). [63 1 The 
magnetic fluctuations in time and space have been described l67l using a localized Heisenberg model with 
competing ferro- and antiferromagnetic interactions between respectively first and second-nearest neigh- 
bors, but to reconcile this model with the partly metallic band structure is a problem. [48 68 69 70, 1\\ 
Another possible solution of the moment problem in SDFT is that moments of the predicted size are 
present, but fluctuate on a time scale faster than what is probed by the experiments. [63 1 In fact, two recent 
studies of realistic, DFT-derived multi-band Hubbard models solved in the dynamical mean-field approxi- 
mation (DMFT) show that the magnetism has two different energy scales. Il72ll73l It is therefore possible 
that the electronic correlations after all do play a role in these multi-band, multi-orbital materials. f74ll75l 

Experiments and calculations have revealed a marked interplay between the details of the band structure 
and the superconducting properties. Most of these observations are empirical and we feel that there is a 
need to explain the origin of such details. In this paper, we therefore attempt to give a self-contained, 
pedagogical description of the paramagnetic and spin-polarized band structures. Specifically, we discuss 
the Fed Asp band-structure topology, causing the pseudogap at d 6 as well as numerous Dirac cones, the 
interlayer hopping in the simple-tetragonal and body-centered-tetragonal structures, the spin-spiral band 
structures, and the band-resolved magnetic energies. In all of this, the covalency between Fe d and As p is 
found to play a crucial role. Applications to superconductivity are beyond the scope of the present paper. 

In Sect. |2]we explain the structure of a single, isolated FeAs layer and use the glide mirror to reduce the 
primitive cell to one FeAs unit and have k running in the large, square Brillouin zone (BZ) known from 
the cuprates. Halving the number of bands will prove important when it comes to understanding the multi- 
orbital band structure. In Sect. [3]we show that this band structure may be generated and understood from 
downfolding, [76 1 of the DFT Hilbert space for LaOFeAs to a basis set consisting of the five Fe d, localized 
Wannier orbitals, or -as we prefer- including explicitely also the three As p orbitals. Even the latter 8x8 
tight-binding (TB) Hamiltonian, h (k) , has long-ranged pp and pd hoppings due to the diffuseness of the 
As p orbitals, and its accurate, analytical matrix elements are so spacious that they will be published at a 
different place. |77| The crucial role of the As p orbitals for the low-energy band structure, the electron 
bands in particular, and the presence of a d 6 pseudogap is emphasized. The different sheets of the FS 
are discussed. In Fig. [2] we show the factorization of the Bloch waves along the lines and points of high 
symmetry in the large BZ. The high symmetry of the single, tetragonal layer allows many bands to cross 
and leads to linear dispersions, and even to Dirac cones. Our understanding of this generic band structure 
of a single layer then allows us to discuss standard DFT calculations for specific materials. This is done in 
Sect.|U where we first see that increasing the As height moves an antibonding p z /d xy level down towards 
the degenerate top of the d xz /d yz hole bands, with which it cannot cross, and thereby causes the inner, 
longitudinal band to develop a linear dispersion. Interlayer hopping is shown to proceed mainly via the 
As p z orbital and to have a strength and (k x , k y ) -dependence which depends on the material family. This 
hopping is strongest for the bet structure where the As atoms in neighboring layers face each other. In st 
SmOFeAs and for k z at the edge of the 3D BZ, the antibonding p z /d xy level reaches the top of the hole 
bands and forms a Dirac cone together with the longitudinal hole band. In LiFeAs and FeTe the Dirac 
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Fig. 2 Upper-right quarter of the large Brillouin zone (BZ) for the glide-mirror space group of the single FeAs 
layer (black), the factorization of the band Hamiltonian (blue), and the LaOFeAs Fermi surface (red). The BZ for 
merely the translational part of the space group has half the area and is folded-in as indicated by the dashed black 
lines. In order to distinguish the corners, M, and edge midpoints, X, of these two square BZs, we use an over- 
bar for the large BZ. Hence r=T, M=X, and X is the common midpoint of the XY and TM-lines. The folding 
causes all three hole pockets to be centered at T and the two electron pockets be centered at M with their axes 
crossed. The blue boxes along the lines of high symmetry contain the orbitals whose Bloch-sums may hybridize 
(belong to the same irreducible representation). At the high-symmetry points, this factorization of the Hamilto- 
nian into diagonal blocks is as follows: F [xy] [XY] [Xz, X] [Yz, Y] [zz, z] , X [xz] [xy, y] [yz, z] [XY, zz, x] , 
and M [XY] [zz] [xy, z] [Xz, X] [Yz, Y] . With often used notations 1431 . the inner and outer sheet of the M-centered 
xz/yz-\ike hole pockets are respectively ot\ and 012 while the X and Y-centered xy/xz and xy/yz-\ike electrons 
sheets are respectively fix and /?2, and the F-centered xy-like hole pocket is 7. 



point is inside the BZ. The interlayer hopping not only causes the As p z -like 2D bands to disperse with k z , 
but also folds the bands into the conventional, small BZ, i.e. it couples h (k) and h (k + 7rx + 7ry). The 
formalism for interlayer hopping is given in Sect. |4.21 and its increasing influence on the band structures 
of BaFe2As2, CaFe2As2, and collapsed CaFe2As2 is shown and explained, for the first time, we believe. 
In CaFe2As2, we find that the nearly linear dispersion of the d xy /p z -like electron band has developed into 
a full Dirac cone. 

The effects of spin polarisation on the generic 2D band structure are discussed in Sect.[5] We consider 
spin spirals which have a translationally invariant magnitude but a spiralling orientation which is given by 
q. Their band Hamiltonian possesses translational symmetry both in configurational and in spin-space, but 
independently of each other as long as spin-orbit coupling is neglected. The spin spiral therefore simply 
couples h (k) to h (k + q) , regardless of whether q is commensurable or not. For h (k) we use the DFT 
pd Hamiltonian derived in Sect.[3] In order to keep the analysis transparent and amenable to generalization, 
we shall treat the exchange coupling using the Stoner model rather than full SDFT. This has the avantage 
that it decouples the band structure and self-consistency problems, so that we can study the band structure 
as a function of the exchange potential, A. In Sect. |5.2| we discuss the bands and FSs for the observed stripe 
order. As long as the moment is a linear function of A, gapping requires matching of d-orbital characters 
as well as FS dimensions (nesting). For larger moments, and ferromagnetic order in the x direction, the 
FS is different and shaped like a two-blade propeller in the k y direction. It is formed by crossing d xy /p z - 
dy Z i and d zz ^j dxy^ bands, which cannot hybridize along the line through the blades and the hub. The 
resulting Dirac cone has been predicted before |78l and also observed. |3ll |79l The interplay between 
pd hybridization and magnetism is discussed using simple, analytical 4x4 models. In Sect. |5.3| we first 
show the static spin-suceptibility, m (A) / A, calculated for stripe and checkerboard orders as functions 
of the electron doping in the rigid-band approximation. The low-moment solution -maybe fortuitously- 
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resembles the behaviour of the observed moment as a function of doping and q. We then discuss the 
electronic origin of the magnetic energies and first show how the magnetic energy may be interpreted 
as the difference between double-counting-corrected magnetic and non-magnetic band-structure energies. 
This directly relates the magnetism to the band structure and we specifically look at the origin of the 
magnetic energy. We find that the magnetic energy gain is caused by the coupling of the paramagnetic d xy 
hole and d xy /p z electron bands, as well as by that of the d xz parts of the two other electron and hole bands. 
The Fermi-surface contributions to the magnetic energy are comparatively small. We can then explain why 
increasing the distance between the As and Fe sheets increases the stripe-ordered moment, and vice versa. 

At the end, we compare our results with those of fully self-consistent SDFT spin-spiral calculations of 
moments and energies as functions of q and doping in the virtual-crystal approximation, for LaOi_ x F x FeAs 
and Bai_ 2 yK2yFe 2 As2. 

2 Structure 

The basic structural unit for the iron-based superconductors is a planar FeAs layer consisting of three 
sheets: (Fig. [T}. In the high-temperature paramagnetic tetragonal phase, the iron atoms form a square 
sublattice (a = 1) with each Fe tetrahedrally coordinated by four As ligands. The latter thus form two 
square lattices above and below the Fe plane at a vertical distance of approximately half the 
a-constant of the Fe sublattice. The Fe and As positions are thus described by respectively : 



where x and y are the orthogonal vectors between the Fe nearest neighbors and n x and n y take all integer 
values, z is perpendicular to x and y, and has the same length. For perfect tetrahedra, 77 = 1, and for 
LaOFeAs, 77 = 0.93. Instead of -q = v / 2cot6'/2 = 2V2zas, it is customary to specify the As-Fe-As 
tetrahedral angle, 8, or the internal parameter, zas- While t are the translations of the Fe sublattice, T = 
nxX + ny Y are those of the As sublattice whose primitive translations are X = y + x and Y = y x. 
The latter are turned by 45° with respect to x and y, and \[2 longer. The translation group of the FeAs 
layer is T and has two FeAs units per cell. These are, however, related by a glide mirror. Rather than 
using the irreducible representations of the 2D translation group, it is therefore simpler to use those of the 
group generated by the primitive Fe-translations, x and y, combined with mirroring in the Fe-plane. These 
glide-mirror operations ("take a step and stand on your head") generate an Abelian group with only one 
FeAs unit per cell and irreducible representations, exp (ik • r) , which are periodic for k in the reciprocal 
lattice, h x 2TTX+h y 2Try, with h x and h y integer. The corresponding Brillouin zone (BZ) shown in Fig. [2]is 
a square, centered at the f -point, k = 0, with corners at the M-points, k = 7ry ± 7rx and — 7ry=F7rx, i.e. at 
±7tX and ±7rY, and edge-centers at the X and Y -points, k = ±7rx and ±7ry. In this paper, we shall use 
this more heavy notation instead of e.g. (tt, tt) for M and (tt, 0) for X as done for cuprates, because for the 
iron superconductors, no consensus exists about whether to use the (x, y) or the (X, Y) coordinate system. 
The overbar is used to designate the high-symmetry points in the 2D reciprocal space for the glide-mirror 
group. In conclusion, use of the glide-mirror group reduces the number of bands by a factor of two, and 
this is important when attempting to understand the intricacies of the band structure. 

In Fig. [3] we sketch the antibonding Bloch sums of the Fe d xy (top) and d xz (bottom) orbitals, and 
realize that with the glide-mirror notation the former has k = and the latter k ■ x = tt. Accordingly, the 
top of the pure Fe d xy band is at T, while the degenerate top of the pure d xz and d yz bands is at M. We 
shall often return to this. (Authors who unfold without reference to the glide-mirror group, may have f and 
M interchanged, with the result that the xy hole pocket and the two xz / yz hole pockets are respectively 
at M and T. In order to avoid this confusion, it is useful to remember that the two xz / yz hole pockets are 
those towards which the electron superellipses at X and Y are pointing). 

The real 3D crystals consist of FeAs layers stacked in the z-direction with other layers intercalated, al- 
though the iron chalcogenides, FeX, have no intercalation. FigQ]specifically shows LaOFeAs, for which all 
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Fig. 3 Sketch of the antibonding Bloch sum of Fe d xy orbitals in the a;y-plane (top) and of the antibonding Bloch sum 
of Fe d xz orbitals in the xz-plane (bottom). A Bloch sum is formed by adding the glide-mirrored orbital multiplied 
by exp (i k • t) , where the glide, t, is a primitive translation, x or y, and the mirror is the Fe plane. The antibonding 
Bloch sum of d xy orbitals has k=0 and that of d xz orbitals has k-x=7r. That the lobes of the real Wannier orbitals 
avoid the As sites (Fig. |6j is indicated by enhancing the countours of the lobes pointing towards the reader. 

our Wannier-orbital (3D) calculations were done, unless otherwise stated. The interlayer coupling is weak 
but not negligible, and it depends on the material. Although the 2D glide-mirror may take the 3D crystal 
into itself, as is the case for LnOFeAs, FeX, and LiFeAs, we do want to use k z to enumerate the states in 
the third direction. For the 3D crystals we shall therefore use the standard 3D translation group according 
to which only the X and Y translations, combined with an out-of plane translation, leave the crystal in- 
variant. The corresponding 2D reciprocal lattice is hx ttX + hy-nY = (hx + hy) 7ry+ (hx — hy) 7rx. 
Hence, the 3D Brillouin zone is as shown by the dashed lines in Fig.[2](for k z =n/2c), with M falling onto 
T and with corners at X and Y, now named M. Interlayer hopping may thus couple the glide-mirror states 
at k with those at k + 7rx + 7ry. This material-dependent coupling will be considered in Section|4]after we 
have explained the generic electronic structure of a single FeAs layer. 

Spin-orbit interaction also invalidates the glide-mirror symmetry, but the splitting of states degenerate 
at k and k + 7rx + 7ry is at most |CFc3d ~ 0.1 eV, and this only occurs if all three xy, yz, and xz states 
happen to be degenerate and purely Fe e?-like. 



3 Paramagnetic 2D band structure 

In this section we shall describe the generic 2D band structure of an isolated FeAs layer. We start by 
observing that the bands are grouped into full and empty, separated by a pseudogap. We then discuss the 
grouping of the bands into Fe 3c? and As 4p, and derive two sets of Wannier orbitals from DFT, one set 
describing merely the five Fe d-like bands and another set describing the eight Fe d- and As p-like bands. 
Armed with those sets, we can return to a detailed description of the low-energy bandstructure, i.e. the 
one which forms the pseudogap at d 6 and the Fermi surface. This is done in subsection !3.3l where we shall 
see that the hybridization between -or covalency of- the As p and the Fe d orbitals is crucial for the band 
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Fig. 4 Band structure of paramagnetic, tetragonal pure LaOFeAs with the experimental structure near the Fermi level 
(= 0) and for k along the high-symmetry lines in the large 2D Brillouin zone (Fig. [2). Band energies are in eV. 
These DFT-GGA bands were calculated with the NMTO method and a basis of Fe d and As p downfolded NMTOs. 
Transformation to real space yields the eight Wannier functions shown in Fig [6] The 2D bands were obtained by 
neglecting the interlayer hoppings and forming glide-mirror Bloch sums of the Wannier orbitals on a single FeAs layer, 
i.e. by appropriately flipping the signs of the intra-layer hopping integrals. The large gaps in the figure are labelled 
by an electronic configuration which corresponds to a p-set and a d-set of Wannier orbitals which span respectively 
the three lowest and the five highest bands. This d-set is illustrated in Fig[5] Upon electron doping in the rigid band 
approximation, the T-centered hole pocket fills once the doping exceeds 0.1 e/Fe, and when it exceeds 0.3e/Fe, also 
the M-centered hole pockets fill. 



topology. Bringing this out clearly, was in fact our original reason for deriving the eight-orbital pd set, 
although the five-orbital d set suffices to describe the low-energy band structure. 

For FeTe and LaOFeAs the formal ionic states are respectively Fe 2+ Te 2_ and La 3+ 2_ Fe 2+ As 3_ . In 
fact, for all parents of the iron-based superconductors, the nominal electronic configuration is ligandp 6 
Fed 6 . The generic 2D band structure is shown in Fig. [4] for energies ranging from 4.5 eV below to 2.5 
eV above the Fermi level and along the high-symmetry lines of the BZ (Fig. 13. In the energy range 
considered, there are eight bands which are seen to separate into three low-energy and five high-energy 
bands. They may be called respectively the ligand p- and iron d-bands, and the corresponding electron 
count is as written on the figure. At p 6 d 6 the two uppermost bands are seen to be detached from the rest, 
except at one (Dirac) point along the XM-line where two bands cross, because their Bloch functions are 
respectively even and odd with respect to reflection in a vertical mirror parallel to XM and containing 
nearest-neighbor As atoms. If the energy of this crossing could be moved up, above the relative band 
maxima at f and M, it would drag the Fermi level along and the material would transform into a zero-gap 
semiconductor. For the iron-based superconductors, however, the Fermi level is merely in a pseudogap 
and the Fermi surface (FS) consists of a T-centered hole pocket, two M-centered hole pockets, and two 
compensating electron pockets centered at respectively X and Y (Fig. |2). 

3.1 Fe d five-orbital Wannier basis 

Characterizing the five upper bands as Fe d is sound, because they can be spanned exactly by five Wannier 
functions 0431 which behave like Fe <i-orbitals. This can be seen in Fig. [5] Our Wannier functions were 
constructed ll76l to have d character on the central Fe site and no d character on any other Fe site. This 
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Fig. 5 The set of five Fe d-like Wannier orbitals (downfolded and orthonormalized NMTOs) which span the five 
LaOFeAs bands extending from -1.8 eV below to 2.2 eV above the Fermi level. Shown are the positive and negative 
contours, x-m (r) = ± |c| , with the former in red and the latter in blue. Orientation and coloring (Fe red, As green, La 
yellow, and O blue) as in Fig.Q] The three orbitals to the left and the two to the right would belong to respectively the 
t% and e representations, had the point symmetry been tetragonal. Now, split into a (d xy ) and t (d xz , d yz ) . Note 
that the t orbitals d xz = (dxz — dy z ) / v2 and d yz = (dxz + dyz) /v2, whose Bloch sums form the proper linear 
combinations for k along FX and XM (Fig.|2]l, are not simply 45° -turned versions of dxz and dyz shown here, in 
particular because the p tails of the latter are on different pairs of arsenics. The p tails are thus always directed towards 
the nearest As neighbors in the same plane, i.e. they are X or Y, and they antibond with the t head. 



makes them localized Wannier orbitals. The five bands of course have characters other than Fe d, and those 
characters are mixed into the Fe d Wannier orbitals. This by-mixing follows the point symmetry in the 
crystal. Specifically, the Fe d xy Wannier orbital has on-site Fe p z character breaking the horizontal-mirror 
symmetry of the pure d xy orbital, as well as strong off-site p z character on all four As neighbors. The sign 
of the As p z character is a/if/bonding to Fe d xy because the As p hybridization pushes the Fe d band up in 
energy. The corresponding nodes between the Fe d and As p tails make neighboring lobes difficult to see 
in the figure. Hence, only the As p z lobes pointing towards the La layers are big. Similarly, the Fe dx z 
Wannier orbital antibonds with px on the two As neighbors in the X direction, and Fe dy z antibonds with 
Py on the two As neighbors in the Y direction. If the Fe-site symmetry had been exactly tetragonal, the 
three above-mentioned Wannier orbitals would have been degenerate and transformed according to the t% 
irreducible representation. However, the non-tetrahedral environment, e.g. flattening of the tetrahedron 
(rj < 1) , increases the energy of the d xy orbital above that of the d orbitals belonging to t, i.e. dx z and 
dy z or, equivalently, d xz and d yz . In LaOFeAs, the energy of d xy is ^0.1 eV above that of dt- The two 
remaining Wannier orbitals, d 3z 2_ 1 =d zz and d y 2_ x 2=dxy, antibond less with As p because their lobes 
point between the arsenics. Fe d zz is seen to antibond with p z on the four As neighbors and Fe dxy 
antibonds with py on the two As neighbors in the X direction, and with px on the two As neighbors in the 
Y direction. In tetrahedal symmetry these two orbitals would transform according to the e representation, 
and that holds quite well also in the real materials where the orbitals are degenerate within a few meV. 
Their energy is ~0.2 eV below that of the dxz an d dy z orbitals. This e-ty. splitting of a central d shell 
in a tetrahedron having p orbitals at its corners is an order of magnitude smaller than the t2 g -e g splitting 
in an octahedron which allows for better alignment of the p and d orbitals. The ^0.2 eV e-t-x splitting 
in LaOFeAs is 20 times smaller than the width of the Fe d-band structure in Fig.@] and does not cause 
separation into two lower e and three higher t-z bands with a pseudogap at c? 4 . Nevertheless, the and e 
orbitals do play quite different roles in forming the band structure near the Fermi level, as we shall see 
later. 

Whereas in cubic perovskites, including the cuprate superconductors, the effective dd hopping in the 
separated ti g and e g bands proceeds almost exclusively through the p tails, which are placed between the 
nearest-neighbor d orbitals, the effective dd hopping in the iron-based superconductors proceeds directly 
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Fig. 6 Same as Fig.[5] but for the set of eight Asp- and Fed-like Wannier functions which span the entire band 
structure shown in Fig.[4] 

between nearest-neighbor d orbitals on the square lattice as well as via the p tails lying above and below 
the plane of the d orbitals. 

Tabulations of the hopping integrals for the d-orbital Hamiltonian may be found in Refs. 1431 and l59l . 

3.2 Fed Asp eight-orbital Wannier basis and its Hamiltonian 

In order to explain the band structure, and in particular the interlayer coupling in Sect. [4] we find it useful to 
exhibit the As p characters explicitly. We therefore choose not to downfold the As p channels, but span the 
entire eight-band structure in Fig.Uby the eight As p and Fe d Wannier orbitals; all other channels remain 
downfolded l76l . These eight orbitals are shown in Fig. [6] Due to the lack of As p tails, the Fe d orbitals 
of this pd set are more localized than those of the d set and the integrals for hopping between them have 
a shorter range. This basis set is also more suited for including the on-site Coulomb correlations. The As 
p orbitals are, on the other hand, quite diffuse and give rise to strong and long-ranged pp and pd hoppings 
inside the layer. For the p z orbital, this is partly due to its La d and O p tails. This situation is very different 
from the one found in the cuprates, where long-ranged pp hopping is blocked by the presence of Cu in the 
same plane. Although the orbitals of the AsFe pd set resemble atomic orbitals more than those of the d set, 
they do tend to avoid the space covered by the other orbitals in the set: The Fe d orbitals avoid the As sites 
and the As p orbitals avoid the Fe sites. This distorts in particular the Fe i 2 orbitals. The on-site energies, 
e Q , and the nearest-neighbor hopping integrals, f Q x ^ y , are given in the table below. Here, all energies are 
in eV and the hopping integral is the matrix element of the Hamiltonian between Wannier orbitals a and j3 
with n x x + n y y being the vector from a to (3, projected onto the Fe plane. All hopping integrals needed 
to obtain converged energy bands together with their analytical expressions will be published in Ref. ifTTl . 
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For LaOFeAs, a=285 pm (and, within a few per cent, the same for the other iron-based superconductors). 
The energies of the p and d orbitals are respectively —1.8 and —0.7 eV. This 1.1 eV pd separation is merely 
a fraction of the 7 eV pd-band width and it therefore seems fair to claim that the band structure is more 
covalent than ionic. Nevertheless, it does split into three lower As p-like and five upper Fe d-like bands 
as noted above. The band structure fattened by the weight of each of the eight Wannier orbitals of the pd 
set is shown in Fig. [7] Here and in the following we write xy for Fe d xy , t for Fe (it, z for Asp z , a.s.o.. 
The strange wiggles of some of the bands may be seen to have strong z character and this tells us that the 
reason for those wiggles is intra-layer hopping via the LaO layers, whose orbitals are downfolded mainly 
into the As z orbital. 
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3.3 2D Bands and Fermi surface 

We now follow the bands around the d 6 pseudogap and begin with the Fermi surface near T. 

The T-centered hole pocket is seen to have xy character and, as sketched at the top of Fig. [3] its Bloch 
function drf-antibonds with all four nearest Fe neighbors. If in this figure we imagine inserting the xy 
orbital of the d set (Fig. q), we realize that the sum of the As p z tails cancel. In fact, none of the other 8 
orbitals in the pd set can mix with the Bloch sum of xy orbitals at f . This we have stated in the caption 
to Fig. |2] together with the selection rules for all other high-symmetry points. The selection rules for the 
high-symmetry lines are given in blue on the figure. For k moving from f towards X, the xy band is 
seen to disperse downwards because in the x-direction, the character of the wavefunction goes from dd 
antibonding to bonding and, at the same time, the band gets repelled by the above-lying x z jy band whose 
xz orbitals point in the direction of the k-vector, i.e. the longitudinal t band. The corresponding inter-band 
matrix element increases linearly with the distance from f , whereby the downwards curvature of the xy 
band is enhanced by about 10%. The resulting hole band mass is about twice that of a free-electron. 

At X, the xy Bloch sum bonds between nearest Fe neighbors in the x direction and antibonds between 
those in the y direction. In addition, weak hybridization with a 2 eV lower-lying As y band provides y 
character to antibond between the xy orbitals in the y direction. This pushes the xy band up at X by 0.2 
eV to —0.4 eV. As k now moves on from X towards M, the Bloch sum of xy orbitals becomes bonding 
between Fe nearest neighbors in the y direction as well, whereby the As z tails (FigO no longer cancel 
and their antibonding contribution increases linearly with k y . This causes the band to disperse strongly 
upwards, to a maximum 1 eV above the Fermi level. As can be seen from Fig. [7] this strong change of 
band character may also be explained as the result of strong pd hybridization and avoided crossing of a 
pure xy band dispersing downwards from X to M and a pure z band dispersing strongly upwards due to 
long-range hopping, partly via La and O (see Fig. [6). Corresponding to the 1.3 eV upwards dispersion of 
the xy/z antibonding band along XM, we see a downwards dispersion the z/xy bonding band. At M, xy 
and z hybridize, but only with each other: the pure xy level is at —2.0 eV, the pure z level at —0.4 eV, and 
the xy-z hybridization is 2 eV, thus pushing the antibonding, predominantly z-like level up to +0.9 eV and 
the bonding, predominantly xy-like level down to —3.4 eV. 

The i-bands which form the M-centered hole pockets exhibit a very similar behavior as the xy band 
when, instead of going along the path T-X-M (or T-Y-M), we go along the path M-X-r for the xz band 
and along M-Y-T for the yz band. The avoided pd crossing is now between a pure xz band dispersing 
downwards from X to f and a pure y band dispersing strongly upwards; these dispersions are strong 
because k changes in the direction of strong hoppings, ddir and ppir, respectively. At X, the xz band is 
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Fig. 7 Band structure from Fig. |4]fattened by the character of each of the eight Wannier orbitals in the pd set (Fig. 
[SJ. A fatness was obtained by perturbing the on-site energy of the orbital in question. That two bands share the same 
fatness means that they hybridize. In order to concentrate the fatness onto as few bands as possible we chose the 
appropriate linear-combination of t orbitals (see Fig. O for each high-symmetry line . The TM-line was chosen as the 
one in the X direction. 



pure and merely 0.1 eV below the Fermi level. At T, the pure xz level is at —1 eV, the pure y level is at 
+0.5 eV, and the hybridization between them is over 2 eV thus pushing the antibonding, predominantly 
y-like level to +2.2 eV above the Fermi level and the bonding, predominanly xz-like level down to —2.7 
eV. Whereas along FX, the xy and xz/y bands hybridize and therefore cannot cross, along XM, they do 
cross because the Bloch sums with k x =ir of Fe xy and Fe xz orbitals are respectively even and odd upon 
reflection in the As-containing vertical mirror perpendicular to the x direction. In other words, they belong 
to different irreducible representations (Fig.|2j. Since outside the XM-line, hybridization between them is 
no longer forbidden, the accidental degeneracy is at a Dirac point in 2D (a line in 3D). To separate the two 
bands would require moving the xy/y above the xz level at X . 

Near X, the xy- and xz-liks bands are close in energy and both have minima. The minimum of the 
xy-like band curves steeply upwards towards M due to strong hybridization with z, and that of the xz-like 
band curves steeply upwards towards f , due to strong hybridization with y, but is flat towards M, which 
is the direction transversal to the dominating ddir hopping. These two bands hybridize weakly with each 
other, except on the XM-line. As a result, they form a lower and an upper band of which the latter cuts the 
Fermi level at an X-centered electron pocket. As long as the Fermi level is well above the band crossing 
along XM, the shape is that of a 4th-order superellipse. This ellipse points towards M where its character 
is mainly xy/z and is flat towards T, where its character is predominantly xz/y. 
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The two bands forming the X-centered electron pocket can be modeled by the Hamiltonian, 

e r „ \ / m z X t~ v xy,xy V xy,xz 

H(k x ,k y )=[ " ( )+( m — w - ^ fc2 y ^ fc2 | (2) 



v xy,xz 



in the basis of two effective (downfolded) Bloch orbitals, xy and xz, and where the origin of (k x , k y ) is 
taken at X. Note that the dispersion of the effective xy band towards M is nearly linear for energies not 
too close to e xy due to the avoided crossing of the pure xy and z-bands and the linear increase of their 
hybridization (see Fig.[7j we shall return to this). In expression (fJJ, e xy =—Q.44eV and e xz =— 0.13 eV 
are the levels at X with respect to the Fermi level, v% y xy =0.53 and v% y xz =—0.29eV-a are band slopes, 
i.e. group velocities, and m% y xy =2.9, 771^^=2.3, and iriddS = —7 are the band masses relative to that of a 
free electron. Negative masses are those of holes. Moreover, r = (ao/a) 2 Ry=0.47eV with ao the Bohr 
radius. The numerical values are for pure LaOFeAs with the experimental structure and were obtained by 
fitting the result of an LAPW calculation near the Fermi level. Like everywhere else in this paper, k is in 
units of the inverse Fe-Fe nearest-neighbor distance, 1/a — 1/285 pm. By symmetry, the xy- and yz-like 
bands give rise to a Y-centered electron pocket which also points towards M, where the character is mainly 
xy/z, and is flat towards f with predominant yzjx character. This is illustrated in Fig. [8]by projection of 
the Fermi surface onto the various orbitals. 

Although none of the sheets of the paramagnetic Fermi-surface have major e character, the two e orbitals 
play a decisive role in the formation of the d 6 pseudogap. The pure zz band is centered at —0.6 eV 
and disperses so little that it lies entirely below the Fermi level. At M, the zz band is pure, and at f, 
hybridization with the z band pushes it up by 0.6 eV, an amount larger than the width of the pure zz band, 
to —0.4 eV. The pure XY band, on the other hand, is broad because the lobes of the nearest-neighbor XY 
orbitals point directly towards each other. This band has its minimum at f , saddlepoint at X, and maximum 
at M. Both extrema, at respectively —1.8 eV and +0.6 eV, are pure. At intermediate energies, the XY 
band is however gapped in large regions centered at X and Y by avoided crossings with the zz band. 

Near the MT-lines, where the XY and zz bands cannot hybridize, the XY band has an avoided crossing 
near M with the upper, transversal hole band, Yz/Y (Xz/X) in the X (Y) direction. The hybridization 
between the transversal hole band and the downwards-dispersing pure XY band vanishes at M, but in- 
creases linearly with the distance from M, and with a slope proportional to the pdir-like hopping integral 
between the XY and Y orbitals. This means that, if at M, the XY level at 0.6 eV could be lowered by 0.4 
eV such as to become degenerate with the degenerate top of the hole bands, then the transverse hole band 
and the XY band would form a Dirac cone. The trace of this cone can still be seen in Fig. [71 in particular 
at the low-energy edge of the fat XY band. If the singly degenerate XY level at M had been below the 
degenerate Yz/Y level, then this lowest level would be the singly-degenerate top of an YY-like hole band 
and the higher-lying, doubly-degenerate level would be the bottom of the transversal Yz/Y electron band. 
The degenerate partner does not hybridize with XY /Y and is therefore independent of the position of the 
X Y-band. Hence, the actual band structure of the LaOFeAs has the XY and Yz/Y levels, at respectively 
0.6 and 0.2 eV, inverted. 

The model Hamiltonian for a Dirac cone is: 




where for the above-mentioned example k is the distance from M. The zero of energy is midway between 
the two M levels, XY and t/p, which are separated by g. In view of the approximately circular shape 
and isotropic XY character of the transversal (outer) M-centered hole sheet seen in Fig.[HJ the isotropic 
2x2 Hamiltonian (0) is a reasonable representation and may be obtained by limiting k to one of the four 
MT directions where only three transversal orbitals can mix (Fig.|2]i, and then downfolding the transversal 
p orbital, i.e. Y if k is along X. Since the unhybridized XY band disperses downwards towards T, its 
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Fig. 8 Fermi surface for LaOFeAs (see also Fig. [2}. The fatness was obtained as in Fig.|7]and thus gives the orbital 
weight times the inverse of the Fermi velocity projected onto the plane. The fatness of the dominating xy, xz, and yz 
orbitals have been reduced by a factor 5. Here, h and e refer to, respectively, hole and electron sheets. 

mass, 77i2, is negative, and since the unhybridized transversal band disperses upwards, its mass, mi, is 
positive. Finally, the coupling vk is proportional to the product of the Yz-Y hybridization at M, which 
is in fact responsible for moving the pure Yz band up by 0.4 eV, and the XY-Y hybridization, which 
is proportional to k. These hybridizations are clearly seen in the pictures of the d Wannier orbitals in 
FigS For this incipient XY—Yz/Y Dirac cone in LaOFeAs, g= 0.4eV, m^— —1-4, mi= 1.4, and v «a 
0.5 eV a = 1.4 cV A = 1.4c/1973, i.e. about thousand times less than the velocity of light, c. Now, for 
k > \g jv = 0.4 — If M| /10, the Hamiltonian © yields the cone: e (k) = ±vk, and for k < \g /v , it 
yields parabolic Yz/Y and XY-\ik& bands gapped by g and with inverse masses given by respectively 
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Later in this paper, we shall meet not only incipient- but real Dirac cones. 

The final gap needed to complete the d 6 gap in the central part of the TXMY square is the one produced 
by the avoided crossings along MT of the downwards-dispersing upper z/xy band with the upwards- 
dispersing upper, longitudinal Xz/ X band. Here again, none of these bands are allowed to hybridize at 
M, and the matrix elements between them increase linearly with the distance from M in the X direction. 
Specifically, the pd matrix elements Xz-z, X-z, and X-xy are all linear in k. Also the zz band at —0.6 eV 
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mixes in, with the weak 2nd-nearest neighbor dd hopping integral between Xz and zz orbitals providing 
the slope of the linear matrix element. At M, the z/xy and Xz/X levels are thus inverted, but by being at 
respectively 0.9 and 0.2 eV, they are too far apart to make the Dirac cone visible in Fig|7] This will however 
change when, in the following section, we consider other materials and include the k z -dispersion. In 
conclusion, the d 6 pseudogap is caused by the XY and uppermost z/xy levels being above the degenerate 
t/p levels at M. Had the opposite been the case, a situation with the two t/p bands entirely above the three 
z/xy and e bands, i.e. that of a d 6 insulator, could be imagined. 

Having sorted out the intricacies of the band structure and thereby understood the subtle origins of the 
X-centered electron pockets and the d 6 pseudogap, we shall finally return to the M-centered hole pockets 
using Figs.[2]|Zl and[8] These hole pockets have fairly complicated shapes and orbital characters. Although 
the t character dominates, p hybridization pushes the top of the band up by 0.4 eV, to 0.2 eV above the 
Fermi level, as has been mentioned before. Departing from M, the two bands split into a steep one with 
relative mass numerically smaller than one and a shallow one with mass numerically larger than one. They 
give rise to respectively the inner and the outer hole pockets. As long as the character of the band is 
predominantly i-like, the steeper, inner band will be the one for which the ddir hopping is along the vector 
distance from M, that is, the longitudinal band. Accordingly, we see in Fig. [7] from M to X, the yz-like 
band stay intact and disperse strongly downwards. From M towards T, we see the Xz-like band disperse 
downwards and stay intact until it suffers an avoided crossing with the zz band. The inner band is in fact 
steeper towards f than towards X and Y, and this is partly because the p hybridization of the longitudinal 
band has a node along Mf , as can be seen for the Xz-like band in Fig. [7] The further reason for the 
small mass of the inner hole pocket is the gapped Dirac cone formed with the 0.7 eV higher-lying z/xy 
band. The outer, transversal hole band has a large mass not given by the weak ddS hopping integral, but 
as discussed above in connection with Eq. (O, by its hybridization proportional to fc, with the g=0.4eV 
higher-lying XY band. Along MT this transversal hole band is mainly xz/x hybridizing proportional to 
fc, not only with XY but also with zz. The latter gives the anisotropy seen in Fig. [8] 

In the next section we shall see how these details are modified by the material-dependent height of As 
above the Fe plane and interlay er coupling. 



4 Influence of As height and interlayer hopping 

Until now we have discussed the generic 2D band structure for an isolated FeAs layer. This band structure 
was obtained by (i) downfolding the proper 3D bands of LaOFeAs with k in the small BZ to a 16 x 16 pd TB 
Hamiltonian, (ii) neglecting the interlayer hoppings and (iii) reducing the resulting Hamiltonian to an 8 x 8 
by transformation to the glide-mirror Bloch representation with k in the large BZ. The understanding of this 
relatively simple, generic, 2D band structure obtained in Sect. l3.3l enables us now to explain the material- 
dependent, complicated, 3D bands obtained by standard DFT calculations in the small BZ. Specifically, 
we shall present and discuss the 3D band structures of simple-tetragonal (st) LaOFeAs and SmOFeAs in 
Fig- E]and in Sect. 14.11 mentioning those of FeTe and LiFeAs en passant, and then in Sect. l4.2l move on to 
the band structures of body-centered tetragonal (bet) BaFe2As2 and CaFe2As2, in the normal as well as the 
collapsed phase. The band structure of bet BaRu2 As2 will finally be mentioned. The interlayer hopping 
is mainly between As z orbitals. In the st LnOFeAs materials this hopping is fairly weak and the material 
dependence of the band structures is caused more by the varying height of As above the Fe plane than by 
interlayer hopping. This we shall see in Sect. 14.11 For st FeTe and LiFeAs, and in particular for the bet 
materials, the interlayer hopping is dominating, and since its effects are non-trivial, we have derived the 
formalism and shall present it in Sect. 14.21 It turns out that the folding of the bands into the small BZ and 
subsequent interlayer hybridization at general k-points cause many bands to have nearly linear dispersions 
and in some cases to form full Dirac cones. 
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Fig. 9 3D band structures of simple tetragonal LaOFeAs (left) and SmOFeAs (middle), as well as body-centered 
tetragonal BaFe2As2 (right). For the two former, the BZ is a rectangular box whose cross-section is the 2D folded-in 
zone shown by the dashed lines in Fig. [2] The Mr-line is in the k z =0 plane, the FZ-line is along the k, -direction, 
and the ZR-line is in the k z =n/c plane. The 2D notation for the projection onto the (k x , fc y )-plane is given on the top. 
The BaFe2As2 band structure is plotted along those same lines, now labelled Xr, TZ,and ZT/2, as may be seen from 
Fig. [10] The computational scheme was GGA-LAPW (|80l). 



4. 1 Simple tetragonal LnOFeAs, FeX, and LiFeAs 

Since in these st crystals, the FeAs and LnRO layers (see Fig. [TJ are simply translated in the z-direction by 
a multiple of c and then stacked on top of each other, the primitive translations in real and reciprocal space 
are respectively 

[ T 2 ) = ( -1 1 | ( y ) and ( £ | = ( -7T 7T ) f y ) . (4) 
V T 3 ) V c A *) V g3 J V * A ') 

The 3D BZ is therefore simply a rectangular box whose cross-section is the 2D zone, folded-in as shown by 
the dashed lines in Fig. [2] The midpoints of the vertical faces, |gi and |g2, are labelled X = (|M, fc z =0) , 
those of the vertical edges, | (gi ± g2) = 7ry and 7rx, are labelled M = (Y, k z =0) and (X, 0), those of 
the horizontal faces, ±^g3, are labelled Z = (f , ±-), those of the horizontal edges, ±i (gi + g3) and 
±i (g2 + g3) , are labelled R =(iM, ±— ), and those of the corners, ^ (±gi ± g2 ± g3) , are labelled 
A=CY,±f)and (X,±f). 

In order to compare with our familiar 2D bands in Fig. |4] we first consider them along XM, where they 
are the same as along YM, and then then translate YM by — g2 to Xf , which is Mr in Fig. [9] Now we can 
easily recognize the M-centered, doubly-degenerate top of the i-like hole bands, the above-lying XY and 
z/xy bands, and the zz band at —0.6 eV. Next, we consider the f X bands in Fig. |4]and translate this line by 
— g2 to MY, which is TM in Fig. [9] This time, we recognize the T-centered xy hole band, the zz/ z band 
at —0.4 eV, and the T-centered bottom of the XY band at — 1.8eV. Near M=X and Y, we also recognize 
the bands responsible for the X -centered electron super-ellipse, along the direction towards f and as well 
as towards M. Finally, we consider the f M bands in Fig. |4]and superpose those translated by gi, to MP, 
onto them. Those bands are symmetric around M/2, and their first half, from f to M/2, is placed from Z to 
R in Fig. [9] Here again, we can easily recognize the bands. After having obtained an understanding of the 
3D band structure of LaOFeAs from merely placing e a (k + g 2 ) on top of e a (k) , we now search -and 
subsequently explain- the effects of interlayer hopping. 

k z is along Mr and ir/c along ZR. Along the vertical path TZ, we see the k z -dispersion at M and f . 
From this, it may be realized that only bands with As z character disperse significantly with k z , i.e. that the 
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interlayer hopping proceeds mostly from As z to As z. The bands seen to disperse in Fig.|9]are the upper 
z/xy band near M and the upper zz/z band near f (Fig.|7]i. This interlayer hopping simply modulates the 
energy of the z orbital, e z (k z ) = t cos ck z . Now we see something very interesting: For k z near tt/c, the 
upper z/xy band has come so close to the top of the M-centered hole pockets that the inner, longitudinal 
band takes the shape of a Dirac cone over an energy region of 0.4 eV around the Fermi level. The inner 
hole cylinder, as well as its radius, thus become warped due to this incipient Dirac cone. This is even more 
pronounced for SmOFeAs because here, the z/xy band lies nearly 0.3 eV lower than in LaOFeAs, in fact 
so low that the 2D band in the k z =n/c plane is nearly a complete Dirac cone at 0.2 eV and with slope 
v ~ 0.3 eV-a (see expression OJ). To bring the Fermi level up to the cusp would, however, require electron 
doping beyond 30%. 

The reason why the z /xy band lies lower in SmOFeAs than in LaOFeAs is that As lies higher above 
Fe (^=0.98) in the former than in the latter compound (77=0. 93) . The z-xy hybridization is therefore 
smaller, and that moves the upper z/xy band down at M. This is clearly seen along all directions in Fig. [9] 
but whereas this flattening of the upper z/xy-like band increases the mass at the X-centered electron 
pocket towards M and makes it more d-like, it decreases the mass of the inner M-centered t/p-like hole 
pocket due to the incipient Dirac cone. Increasing 77, generally decreases the pd hybridization, whereby 
pd antibonding levels move down in energy with respect to those of pure d character and become more 
d-like. Important effects of this are the lowering of the top of the t/p hole band at M with respect to that 
of the xy hole band at T and the lowering of the bottom of the xy-like electron band with respect to the 
pure xz level at X . These changes are clearly seen in Fig. [9] With the Fermi level readjusted, the size of 
the f-centered xy hole sheet is increased for SmOFeAs and is now similar to that of the outer M-centered 
hole sheet. Finally, we may note that the decreased t/p hybridization at M decreases the coupling linear in 
k to the XY band, so that this band becomes less steep in the Sm than in the La compound. 

Of all known iron-based superconductors, SmOFeAs has the highest T c max (55 K) and the most regular 
FeAs4 tetrahedron, i.e. its 77 is closest to 1 . For nearly all Fe-based superconductors, T c max versus 77 seems 
to follow a parabolic curve, |[T3l a correlation which has been extensively studied, but is not understood. 
For LaOFeAs, T c max =27K. 

Also LiFeAs and the iron chalcogenides, FeX, have the st structure and calculations l23l [34l [35ll yield: 
77=1.12 for LiFeAs, while for X = S, Se, Te: ?/=0.87, 0.97, 1.16, respectively E3- In LiFeAs and FeTe, 
the upper z /xy band thus sits considerably lower in energy. Moreover, since in LiFeAs the perpendicular 
As z hopping is enhanced by hopping via Li s, and since the perpendicular Te hp z hopping is stronger than 
the As 4p z hopping in LnOFeAs, the z/xy band disperses 3 and 4 times more along TZ in respectively 
LiFeAs and FeTe, than in LnOFeAs. As a consequence, the z/xy -like band crosses the degenerate t/p 
band already when k z ^n/2c, and here, it forms a Dirac cone with the inner, longitudinal t/p band, at 0.1 
eV above the Fermi level in LiFeAs and at 0.2 eV in FeTe. For 7r/2c < k z < 37r/2c, the band which at 
M has longitudinal t/p character disperses upwards and the other band, which at M has z/xy character, 
downwards. Accordingly, the inner hole sheet of the Fermi surface is not a cylinder, but extends merely a 
bit further than from — it /2c to 7r/2c where the z/xy band along TZ dips below the Fermi level. The mass 
of this sheet vanishes when k z is at the Dirac value, ~7r/2c. Here, the slope of the cone in the (k x , k y )- 
plane is v ~0.5 eV-a. The T-centered xy pocket is a straight cylinder, whose cross-section in FeTe has 
about the same size as that of the outer M-centered t/p hole sheet at k z =0, i.e. like in SmOFeAs, and in 
LiFeAs is even a bit larger. LiFeAs is a non-magnetic superconductor with T c =18 K. 

4.2 Body-centered tetragonal BaFe2As2, CaFe2As2, and BaRu2As2 

In the body-centered tetragonal (bet) structure, the FeAs layers are translated by x before they are stacked 
on top of each other. This means that the As atoms of adjacent layers are directly on top of each other. 
Moreover, the interlayer As-As distance, d =379 pm is about the same as the intralayer As-As distances, 
pm and \[2ar\. It is therefore conceivable that the interlayer hopping vertically from As z to As 
z (ppa) is very strong. This is in fact the reason for the 2 eV dispersion seen along TZ in the band structure 
of BaFe2As2 on the right-hand side of Fig. [9] Since 77=0. 97 for BaFe2As2, the position of this z-like band 
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is not as low as in FeTe, but more like in SmOFeAs. At fc 2 «37r/4c, the band crosses the degenerate t/p 
band and forms a Dirac cone with its longitudinal branch in the (k x , fcj,)-plane, as we shall see explicitly 
later. Note that the longitudinal branch disperses downwards from V towards X in the k z =0 plane, but 
upwards from Z in the k z =7r / c plane. 

Ba is intercalated in the holes between the neighboring As sheets and thus has 8 nearest As neighbors. 
Also Ba orbitals can be vehicles for interlayer coupling and, in fact, a Ba hd xz / yz band lying above the 
frame of Fig.|9]repels the top of the doubly degenerate t/p band near T with the result that there, the latter 
is only slightly above the top of the dispersionless f xy band, whereas at Z, it is 0.2 eV above. Clearly 
visible in the figure is also a Ba 5d xy band starting at 1 .0 e V at X and then dispersing downwards towards 
r, which is reached at 0.3 eV after an avoided crossing with the y/xz band decreasing from its maximum 
at r. From T towards Z, the 5d xy band then disperses upwards to 1.0 eV and, from there, continues in the 
k z =ir/c plane towards T, but soon suffers an avoided crossing with the hybridized z/xy-longitudinal-i/p 
band. 

But before we continue our discussion of the BaFe2As2 bands we need to write down a formalism for 
the interlayer coupling which is strong -and poorly understood- in the bet structure. 

We start from the 2D Bloch waves, \r;a, k) , of a single FeAs layer with a labelling the state (e.g. the 
band) and k the irreducible representation of the glide-mirror group. These 2D Bloch waves are expressed 
as linear combinations of localized Wannier orbitals. For the 3D crystal, we now use its out-of-plane 
translations, 113T3, to stack the 2D Bloch waves in the 3rd direction and form the corresponding Bloch 
sums: 

l r ;a> k )=y] |r-n 3 T 3 ;a,k) exp(m 3 T 3 -k), (5) 

— '713=0, ±1,... 1 

which we shall then use as basis functions. Here and in the remainder of this chapter, an overbar is placed 
on the 2D Bloch vector in order to distinguish it from the 3D one, k = k+fc z z. 

Since in the bet structure, the As atoms in a top sheet are vertically below those in the bottom sheet of 
the layer above, the corresponding vertical interlayer hopping via As z is particularly simple and strong, so 
we first specialize to this case. The bet primitive translations in real and reciprocal space are respectively: 
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and g 2 = -TT 7T -f y , (6) 

V S3 J { * A •) 

where c is the distance between the FeAs layers in units of the nearest-neighbor Fe-Fe distance, i.e. half 
the bet c-lattice constant. The bet Brillouin zone and the stacking between neighboring Wigner-Seitz cells 
(BZs) of the reciprocal g-lattice is shown in Fig. [10] Using that T3 ■ k = k x + ck z , it is now a simple 
matter to form the 3D Bloch sums ([5]). For two states a and f3 with the same k, the interlayer coupling 
caused by the vertical z-z hopping is easily found as: 



H 



P,k) = t 1 - c* (k) c 2 ,£ (k) coscfc z , (7) 



inter 



where t 1 - is the ppa hopping integral (< 0) between an As z orbital in the top sheet to the As z orbital 
vertically above, in the bottom sheet of the next layer. c Zi/ g (k) is the eigenvector coefficient to the As 
z orbital in the 2D /3-state. Note that we have not missed a factor 2 in (|7), because only one lobe of 
the z orbital is used for interlayer coupling. Now, from the mere knowledge of a k-function in a single 
-bottom or top- As sheet, the k and k + g 1 (or k + g 2 ) translational states are indistinguishable. Their 
difference is that they have opposite parity upon the glide-mirror interchanging the top and bottom sheets. 
Interlayer hopping can therefore mix states with k and k + g 1 . For the corresponding interlayer coupling, 
qj(f -diagonal in the 2D Bloch vector, we then find : 



H 



■Si) = -^c* (k) c Zj/3 (k + gx) sincfc,,, (8) 

';er 
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Fig. 10 Central bet Brillouin zone and the one translated by the reciprocal-lattice vector gi = 7rx + ny + nz/c, 
see Eq. J6). V is the center of the zone, Z = ±g3/2 = ±7rz/c are the centers of the 2 horizontal faces, N are those 
of the 8 large slanting faces, e.g. gi/2, X are the centers of the 4 vertical faces, e.g. (gi — g2 + g3) /2 = 7rx and 
(gi + g2) /2 = ny, and P are the 8 corners between the vertical X neighbors, e.g. 7rx + 7rz/2c. Note that Z in the 
zone translated by gi is gi — g3/2 = 7rx + ny. 



where the different parities of the k and k+gi states causes the sin ck z -dispersion. Finally, before coupling 
the k + gi state (see Fig. |T0T > to that with k, the former must be brought back to the central zone, and that 
requires shifting k z in (0 back by ir/c. As a consequence, 



,k + g x 



H 



/S.k + gA =-t- L < Q (k + g 1 )c Zi/j (k + g 1 )coscfc z . (9) 

/ inter 



In order to get a first feeling for this formalism, let us assume that we have nothing, but interlayer 
hopping. That is, we have pure z states which only couple between -but not inside- the layers. The 
Hamiltonian for this problem with two z-orbitals per cell is: 



H(k) = 



t cos ck z —it sin ck z 
it^-smck-, — t^~ cos ck? 



Diagonalization yields two dispersionless bands with energy ±i x , and this is because this system without 
intra-layer coupling is merely an assembly of As2 dimers (dangling bonds). The same kind of thing hap- 
pens at the non-horizontal boundaries of the bet BZ, where the states with energies e a (k) and e a (k + g-J 
are degenerate, because if there are no further degeneracies and if t is so small that we only need to con- 
sider those two states, their energies simply split by it" 1 |c 2 Q (k) | . We thus see, that neglecting inter- 
layer coupling does not simply correspond to taking k z =ir/2c; this merely makes the diagonal couplings 
vanish. 

In simple tetragonal FeX, the interlayer coupling proceeds mainly from an As z orbital to its 4 nearest 
As z orbitals in the next layer, with a hopping integral i z . From this follows that the diagonal and off- 
diagonal interlayer couplings are given by respectively 



a, k 



H 



/3, k) = At c* z a (k) c Zy p (k) (cos k x + cos k y ) cos ck z 



inter 



H 



/ 8,k + g 1 )_ = — Ait c* a (k) Cz,/3 (k + g-J (cosk x + cosk y ) sincfc z . (10) 



Since the simple tetragonal reciprocal lattice vectors gi and g2 in Eq. (HJi have no k z component, there is 
no k z -translation leading to an Eq. (9). Instead, the prefactor cos k x + cos k y provides the sign-change for 
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k going to the next BZ. In addition, the prefactor makes the interlayer coupling vanish on the vertical faces 
of the 3D st BZ, that is on (XY, k z ) . This is different from the bet case. Finally, for the system with only 
z orbitals and no intra-layer coupling the st fomalism yields two dispersionless bands with energies ±4£ z 
corresponding to isolated AS-AS4 molecules with no coupling between the 4 atoms in the same plane. 

Having deepened our understanding of the interlayer coupling via the As z orbitals, we can now return 
to our description of the bet band structures for which this interlayer hopping is particularly strong and - 
most noticeably- gives rise to the 2 eV dispersion of the z-like band seen along TZ for BaFe2 As2 in Fig. [9] 
This dispersion is seen to be five times larger than in the LnOFeAs compounds, and it turns out that this is 
not even the entire interlayer dispersion, \c zz / xy (M) | : 

From our previous discussion of the 2D bands, we may recall that, at M, the z orbital can only mix 
with xy and the level of interest is the antibonding z/xy level which in LaOFeAs is at 0.9 eV and has 
about 70% z character, specifically, c zz / xy (M) = — 0.83. At f , z can only mix with zz and there again, 
the level of interest is the antibonding zz/z level, which in LaOFeAs is at —0.4 eV and 85% d-like, 
c z,zz/z (?) = — 0.39. If we now just couple those two LaOFeAs antibonding bands with the bet interlayer 
coupling given by Eqs. ©, (8), and (0, and adjust the one parameter t to the BaFe2As2 TZ band, thus 
yielding t ~ —2 eV, we get the two bands shown in Fig.QT] The good agreement of the upper band with 
the TZ band in BaFe2As2, dispersing from 2.1 to 0.1 eV in Fig. [9] hints that this band does result from 
an avoided crossing of the downwards-dispersing Mz/xy band and the upwards-dispersing T z/zz band, 
such that the 0. 1 eV state at Z is not M zj xy, but f zz/z. The M z/xy state at Z must then be the top of the 
lower band which is seen to have energy —0.9 eV in BaFe2 As2 (this includes a push-down by a high-lying 
Ba hd zz band). The bottom of the lower band is then the f z/zz state at T, which in BaFe2As2 is seen 
accidentally also to have energy —0.9 eV. 

After this estimate, let us briefly recall the proper way of including the interlayer z-z coupling along 
TZ. First of all, from the caption to Fig. [2] we learn that z can only mix with xy at M and with zz at T. 
Secondly, from the 2D bands in Fig. [7] we see that at M, as well as at T, the z-like levels are separated 
by as much as 4eV. This is due to strong pd hybridization; the z, xy matrix element at M is 2 eV and the 
z, zz element at f is 1.2 eV. For an interlayer hopping, t , as large as 2 eV, we now ought to solve a 4 x 4 
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T k =n/c M M 



r k z =m/c Mm y r rc/2c M M 




3k/4c M 



Fig. 12 3D band structure of bet BaFe2 As2 , fattened by various partial-wave characters and calculated with the 
GGA-LAPW method ([80)). The k= (k,fc z ) paths chosen are: (f XMf, 0) = (MYfM, tt/c) , in the two top left 
panels, (f XMf , n/2c) = (MYfM, n/2c) in the two top right panels, and (fXMf,7r/4c) = (MYfM, 3tt/4c) 
in the four bottom panels. Note that these paths extend over the two BZs shown in Fig. [10] The fatness of As z has 
been enhanced by a factor 5 compared with those of Fe d in order to account approximately for the fact that the As z 
Wannier orbital has most of its charge density outside the atomic sphere used in the LAPW method. 



eigenvalue problem (that the d states at k and k + g are different is irrelevant). However, for producing 
Fig-CU we got away with neglecting the pd-bonding levels. On the other hand, t ~ —2 eV was obtained 
by fitting to bands which include the Ba 5d zz hybridization and this leads to an overestimation of t- 1 , as 
we shall see later. 

We now discuss further aspects of the bet band structure computed for BaFe2As2. The first two top 
panels of Fig. [12] show respectively the As z and Fe zz projected bands for k z =0 and along the path f X 
Mr, familiar from Fig.|7]but labelled TXZT in the bet reciprocal lattice (Fig. [Tot . This path takes us from 
the center, T, of the central BZ to the center, X, of a vertical BZ face, from there to the center, Z, of the 
bottom face of the neighboring BZ, and finally back to the origin, V. The piece outside the central BZ may 
of course be translated back to the bottom face of the central zone by — gi , or to the top face by g3 — gi . As 
a result, not only the eight TXZr= (f XMf , 0) -bands, but also the eight TXZr= (MYfM, tt/c) -bands 
are obtained in such a standard calculation. For k z =Q and tt/c, the k-states are pure because smck z =0, 
and the Fe zz projection confirms that the 0.1 eV state at Z has zz character. 

The other band with zz character, i.e. the M zz band, is seen to have energy —0.6 eV at V and — 1.0 eV 
at Z. For reasons of symmetry, this band can have no As z character, and nevertheless disperses with k z (as 
seen directly along TZ in Fig.|9]l. This is due to repulsion from the Ba 5d xy band. 
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When, in Fig. [9] we observed that the longitudinal branch of the doubly-degenerate t/p band curves 
downwards at V in the fc 2 =0 plane, but upwards at Z in the k z ~n/c plane, both in the central BZ (Fig.llOb. 
this was seen as a consequence of the Dirac cone in the fc 2 =37r/4c plane at the crossing of the t/p band 
with the upper z-like band between T and Z. In the meantime, we have learned that the upper z-like band 
cannot hybridize with the t/p band in the k z =n/c plane where s'mck z =Q, because the former is a T-state 
and the latter an M-state at Z. The upwards curvature at Z is therefore due to repulsion from the lower 
z-like band, the M z/xy band near —0.9 eV. This repulsion is substantial and causes another, but merely 
incipient Dirac point at —0.4 eV at Z. However, as we now depart from the k z =nir / c plane, also the F z/zz 
band hybridizes with the longitudinal t/p band -proportional to sin ck z and to the horizontal distance from 
M- and this causes a Dirac cone to be formed at the crossing of the upper z-like and the degenerate t/p 
bands at k z =37r / 4c and 0.25 e V. At the bottom four panels of Fig. [T2]we therefore show the band structure 
in the k z = ir/Ac + nir/c and k z = Sir /4c + nir/c planes extending over the two BZs, projected onto the 
As z and relevant Fe d partial waves. The Dirac cone at (M,7r/4c) in the BZ at gi is seen to have z, zz, 
longitudinal t/p, and xy characters in agreement with what was said above, and to have a low-energy slope 
and a steeper high-energy slope. The low-energy slope v = 0.4 eV a is of course due to the hybridization 
of the t/p band with the crossing, upper z-like band and the high-energy slope v = 0.8eV-a is due to 
hybridization with the lower z-like band, which is around —0.9 eV. (The hole band seen in Fig. [l2]to have 
the strongest xy character is irrelevant for the Dirac cone because it is the f -centered hole band which, due 
to lack of z character, cannot mix with k + gi states). As k z is now increased above 7r/4c, e.g. tofc z =7r/2c 
as shown in the last two top panels of Fig. [12] the upper z-like band moves above the t/p band at M to 

1 .0 eV, whereby the longitudinal branch of the latter curves downwards, like the transversal branch. For 
k z ='Sir j '4c, the behaviour can be seen near M in the four bottom panels: The upper z-like band is now at 

2.1 eV and has no z/zz but only z/xy character. The repulsion from this band steepens the longitudinal, 
inner t/p hole band. 

The outer, transversal t/p hole band attains its hole character mainly from repulsion by the XY band, as 
was explained in connection with Eq. (0 and seen in Fig. |7](but excluded in Fig.[T2l. Finally, for k z ~n/c, 
the behaviour can be seen near M in first panel, which in fact shows what we have already observed without 
orbital projections on the right-hand side of Fig. [9] 

The hole part of the Fermi surface thus has a strongly warped, cylindrical sheet, centered around the 
vertical TZ line (see Fig. [Tot . In most of the zone, this sheet has longitudinal Mt/p and some M z/xy 
character, i.e. it is the small-mass, inner hole sheet. Going from T towards Z, this sheet narrows down to a 
neck for the Dirac value, k z =kp z , where the T z/zz character starts to dominate. Finally, close to Z, the 
character becomes purely T z/zz and the cylinder bulges out. Unlike in FeTe and in LiFeAs, this sheet 
remains a cylinder because the upper z-like band remains above the Fermi level for all k z . Concentric with 
this longitudinal t/p hole cylinder is the transversal one, whose mass is dominated by its XY character. 
That cylinder has little warping and lies outside the longitudinal cylinder, except near Z where the latter 
bulges out. The third hole sheet has nearly pure T xy character and is a straight cylinder. Also this cylinder 
is centered along TZ, but it does not hybridize with the two M cylinders and it is as narrow as the Dirac 
neck of the longitudinal t /p cylinder. 

Next, we turn to the electron sheets, which in 2D are the X and Y-centered super-ellipses pointing 
towards M (Fig.|2]i. They are formed by the xy-Yike band together with the transverse xz-like band for 
the X sheet and with the transverse y z-like band for the Y sheet (see Eq. (O). Figs.|7]and[8]show that 
the only part which has As z character and may therefore disperse with k z , is the one pointing towards 
M. The corresponding band is the z/xy band which we studied above and which at M was found to 
disperse by 3 eV, from 2.1 eV (T) to —0.9 eV (Z). Being centered at respectively X and Y, the electron 
cylinders are however far away from M, and since the xy band can have no z character at X and Y, it 
hardly disperses with k z there. But the z character increases linearly with the distance k y from X -and 
with the distance k x from Y- towards M, and so does the upwards k z -dispersion. As a consequence, the 
latter is strongest towards T= (M, 7r/c) and weakest towards Z = (M,0), where the z/xy band eventually 
bends over and becomes part of the longitudinal t/p band dispersing upwards from Z (see bottom panels 
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of Fig.[T2t. This diagonal interlayer coupling (|7]i thus modulates the long axis, 2/cf, of the super-ellisoidal 
cross section such that it becomes minimal towards T and maximal towards Z (see Fig. [Tot . Specifically, 
for the long axis of the Y cylinder: 2kp x (k z ) w 2kp+8kF cos ck z , and the same for the long axis of the X 
cylinder, 2kp y (k z ) . Taking then the coupling of the X and Y cylinders into account, we first translate the 
X cylinder to the Y site and upwards by n/c, and then couple the two cylinders by the matrix elements (0. 
The coupling has no effect in the k z =nir/c planes containing the T, X, and Z points, but in the k z =n/2c 
planes containing the P and N points, the two cross sections mix to become identical around P, as can be 
seen in the last two top panels of Fig. [12] As a result, the double cylinder twists and follows the shape of 
the string of X-centered rhombic BZ faces (Fig.|2]i, i.e. it stretches out towards the zz /z-bulge around Z of 
the longitudinal hole sheet. Finally we should mention that the other part of the electron double cylinder is 
made up of the Xxz/y and Y yz/x bands which have no z character and no k z -dispersion. 

In the two last two top panels of Fig.[T2] we observe a Dirac point at P and — 1.4 eV. Its upper cone 
in this k z =n/2c plane stays intact over an energy range of nearly 1.5 eV and over a distance of almost 
7r (v = 0.5eV-a) whereafter it develops into the 4 neighboring maxima of the mixed M t /p-longitudinal 
and T zz/z hole band. The lower cone cone extends merely over 0.2 eV. In addition, the Dirac point 
has a second, upper cone which slopes by as much as v = 1.3eV-a and extends several eV above the 
Fermi level, but is truncated slightly below. This is the xy/z electron band. From the cross sections of 
the band structure with the planes shifted by multiples of 7r/4c in Fig. [12] one can see that the two upper 
Dirac cones at the P points are fairly 3D. This seems to differ from the previously discussed Dirac point 
at 0.25 eV and (M, fcn z ) whose cones merely extend in a particular k z =kj^ z 2D plane. In that case, the 
mechanism is that k z tunes the relative position of two bands, which have different symmetries at a 2D 
high-symmetry point, kn, and a hybridization increasing linearly with the distance from that point, to 
be degenerate at k = (kn, k^ z ) ■ Referring to expression (f3): k z tunes the gap, which vanishes at the 
Dirac point, g (k z —k-D Z ) — 0. For the low-energy z/zz-t/p cone at 0.25 eV and (M, 37r/4c) , g (k z ) is 
the k z dispersion of the upper Mz/xy-Tzz/z band seen in Fig. [TT] which is larger than that of the cone, 
v = 0.4eV-a. For the Dirac cones at P, the mechanism is actually the same, but the k z dispersions of 
the two relevant z-like levels are merely from —1.8 eV and — 1.4 eV at P to — 1.6eV and — 1.2eV at X 
(Fig.[T2l. and this amounts to less than the Dirac slopes, v = 0.5 and 1.3 eV-a. 

The reason for the small k z -dispersion along XPX, compared with that along TZ, is simply that XPX is 
at the zone boundary where interlayer coupling is between the degenerate X and Y states. Hence, to first 
order in t- 1 , the two degenerate a levels split by |c ZiQ (X) | , which is independent of k z . But if t x 
were as large as 2 eV, first order might not suffice. Therefore, once again, we first consult the caption to 
Fig.[2]to learn that at X , z may only mix with one other state, yz, which is the bottom of the M longitudinal 
t/p band. The interlayer coupling of X with Y therefore requires merely the solution of a 4 x 4 matrix, 
which after exact Lowdin downfolding of the d block reduces to: 

e ^ fc ') =e '^ + e( x,S (x ! <i (x) ±f± - (n) 

These two second-order equations for e (X, k z ) can be solved exactly and give no k z -dependence. The 
above-mentioned 0.2 eV k z -dependence seen in Fig. [T2] must therefore be due to interlayer hopping via 
orbitals other than As z, but this is negligible. Secondly, we confirm from Fig. [7] that at X , there is 
essentially only one level with z character (because t 2 J)Z (X) = 0.1 eV). This level is at —2.1 eV and is 
non-bonding between nearest neighbors separated by x±?/z or y ± r/z. The yz level is at —1.3 eV and is 
the one seen in Fig.[T2]for BaFe 2 As 2 to be at — 1.4eV at P and — 1.2 eV at X. So if BaFe 2 As 2 were merely 
2D LaOFeAs with added bet interlayer coupling, the shift from —2.1 to — 1.6eV should be ~ \t~ L \ , but 
that is too inaccurate. In order to find t , we therefore seek the lower z-level along XPX in BaFe 2 As 2 and 
find that it is dispersionless and lies 2.0 eV below the upper z level (and below the frame of the figure). 
Hence t = —1.2 eV, assuming a 20% bymixing of yz character to the upper z-like level due to the level 
separation of merely 0.4 eV in BaFe 2 As 2 . 
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Fig. 13 3D band structure of bet CaFe2As2, fattened by the As z character and along the paths (rXMr, 0) = 
(MYfM, 7r/c) and (fXMf,7r/2c) = (MYfM, 7r/2c) . First two panels: normal pressure. Last two panels: 
Collapsed phase at 0.48 GPa.(82) Otherwise as Fig. [12] 

We remark that although the bands do not disperse along XPX, the wave-function characters of course 
do, e.g. in the k z =mr/c planes, the upper and lower z-like levels have respectively purely X z/yz and Y 
z/xz characters, while in the k z =ir/2c plane, they are completely mixed. Secondly, we remind that the z 
states along XPX are intra-layer non-bonding and interlayer ppa bonding and antibonding. 

Going away from the XPX line, the upper z-like band hybridizes linearly with the nearby zz/z and 
z/xy bands and thereby form the P-centered Dirac cones discussed above. That the cones are centered at 
the k z =n/2c plane, i.e. the one containing the high-symmetry points N and P (see Fig. [Tot, is due to the 
fact that, in this plane, the bet bands are periodic for k in the small BZ, whereby the bands along TX equal 
those along MY, which by tetragonality equal those along MX . This higher symmetry is clearly seen in 
the last two top panels of Fig 0] 

In CaFe2As2, 77=1.04 so that the intra-layer z/xy and z/zz hybridizations are smaller than in Ba, 
whereby the centers of the antibonding z jxy and z/zz bands lie lower. However, of greater importance 
is that the smaller size of the Ca ion makes the As-As interlayer distance 70 pm -i.e. nearly 20%- shorter 
than in in BaFe2As2. This substantially increases the interlayer hopping t^ . As a result, the splitting of the 
z-band at X, ~ 2t , which was 2 eV for Ba, is 3.2 eV for Ca. As seen in the first two panels of Fig. [13] 
this splitting is from —3.9 to —0.7 eV. Hence, the upper z-level is above the X yz and zz/z levels at X 
and merely 0.1 eV below the xy-like level at —0.6 eV. The consequence is that in Ca, the xy-like electron 
band forms a complete, upper Dirac cone in the k z =ir/2c plane. This cone slopes by about 3.5 eV over 
the distance tt (v = 1.1 eV • a) and is thereby twice as steep as the Yz/XY cone considered after Eq. ([3J. 
Contrary to the case in BaFe2As2, the xy-like band now forms the inner electron cylinder. 

The second consequence of the upper z-level at X lying as high as — 0.7eV, is that the f z/zz hole 
band lies higher than in Ba. At Z, this amounts to 0.4 eV and places the f z/zz level 0.2 eV above the 
degenerate Mt/p level. Going out in top face of the central BZ fromZ, this f z/zz hole band now crosses 
the longitudinal t/p electron band which, like in BaFe2As2, curves upwards due to repulsion from the 
Mz/xy band (hybridized with Ca4d 2Z ). This high-lying T z/zz hole band crosses the Fermi level far 
outside the transverse Mt/p hole band. For k z decreasing below tt/c, the crossing between the T z/zz 
hole band and the M z/xy-hybridized longitudinal Mt/p electron band gap proportional to sin ck z . The 
resulting lowest band is thus shaped like a volcano with a wide foot of f z/zz character and a caldera of M 
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t/pz/xy character around Z and 0.5 eV above the Fermi level. As k z decreases towards 37r/4c, this gap 
increases so much that the rim is washed out and the caldera develops into a flat hilltop. Eventually, the 
characters of the z-like bands (Fig. [TTT i gapped around the degenerate M t/p band along ZT change back 
to normal order with the upper band being Mz/xy like and the lower band T z/zz like. As a consequence, 
the flat hill continuously transforms into the inner, longitudinal Mt/p sheet. This can be seen for k z ~n/2c 
in the second panel of Fig. Qj] Since the two z-like bands are gapped around the Mt/p band along ZT, 
there is no Dirac point along TZ. The corresponding sheet of Fermi surface is thus a warped, TZ-centered 
cylinder with a very broad, f z/zz-like base near Z and a long narrow piece around T. Outside of this, 
except near Z, lies the transversal Mt/p hole cylinder. The top of the T xy hole band is slightly above that 
of the Mt/p band and its straight cylindrical FS sheet lies outside the transversal Mt/p hole cylinder. 

Pure CaFe2As2 becomes superconducting with T cmax ~12K without doping but by the application of 
hydrostatic pressure in the range 2 — 9kbar. l83l At 5.5 kbar there is a first-order phase transition into a 
collapsed bet non-magnetic and possibly superconducting phase l84l with rj marginally smaller and with 
the interlayer As-As distance decreased by an additional 27 pm. In this collapsed phase, whose bands are 
shown in the two last panels of Fig. Qj] the interlayer hopping is increased so much that the splitting of the 
z-band at X, ~ 2t J - , is now 4.5 e V. As a result, the upper level is 0.2 eV above the Fermi energy. This, first 
of all means that the electron cylinder has lost one sheet, essentially the xy / z sheet, so that there is no Dirac 
cone at P. On the other hand, at Z, the f z/zz level is now 1 eV above the Fermi level and the Mz/xy level 
is 0.1 eV below the doubly-degenerate Mt/p level. The latter creates a prounced, slightly gapped Dirac 
cone with v = 0.8 eV-a. Moreover, the doubly degenerate t/p level is essentially at the Fermi level. As k z 
decreases below tt/c, the T z/zz hole band and the upper Dirac cone, which has mixed longitudinal t/p 
and z/xy character, gap at their crossing, which is at 0.7 eV. This volcano thus has a cone-shaped caldera. 
As k z decreases towards 7r/2c, the rim and the Dirac caldera are flattened away and this flat hilltop sinks 
below the Fermi level. The holes are thus in a large f z/zz M xy/z like sheet shaped as a disc centered 
at Z. At this center, there may be a non-occupied pin-hole with Dirac character. There are no f xy holes 
because that band is slightly below ep. As mentioned above, the electrons are in an XP-centered cylinder 
of xz/yz character. This band structure is thus very different from the standard one, but quite interesting. 

The band structure of BaRu2As2 ll85l is similar to that of (non-collapsed) CaFe2As2 shown in the first 
two panels of FigQj] including a Dirac point at P. But it differs in two respects: The f xy hole band is 
entirely below the Fermi level and the doubly degenerate top of the hole bands disperses like in BaFe2As2 
due to hybridization with Ba 5d xz / yz near V. This causes the top of the t/p bands to sink below the Fermi 
level near T and the corresponding inner and outer Fermi-surface sheets to truncate. 

5 2D Spin-spiral band structure 

At low temperature and normal pressure, the parent compounds of the Fe-based superconductors (except 
LiFeAs) become orthorhombic, antiferromagnetic metals. Superconductivity seems to appear, once these 
spin and charge orders are suppressed, e.g. by doping (electron or hole) or pressure. This superconductivity 
is presumably mediated by spin-fluctuations. 

In this section we shall study the interplay between the band structure and the spin order. Ab-initio 
calculations employing spin-density-functional theory (SDFT) tend to yield the proper spin order, which 
is striped with spins on the iron rows along x aligned and along y alternating. The moments, albeit still 
considerably smaller than the saturation magnetization of 4 fi s/Fe, are much larger (^2 /Fc) than those 
obtained by neutron scattering or muon spin rotation; the latter are ~^0.4 [is /Fc for LaOFeAs fTZI and 
twice as large for BaFe2As2. E] Even worse, only with the calculated large moments do the spin-density- 
functional calculations yield the correct structure (77=0. 93 and 0.5% contraction in the ferromagnetic di- 
rection for LaOFeAs); the structures calculated without allowing for spin-polarization differ much more 
from the observed ones than is normal for density-functional calculations (?y=0.81 and no orthorhombicity 
for LaOFeAs). Il38ll63l It thus seems that the large moments exist, but fluctuate on a time scale shorter than 
what can be resolved with neutrons or muons f72l . Below, we shall first study how the spin-polarization 
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modifies the band structures discussed in the previous section which were calculated for the experimental 
structures. Thereafter we shall consider the energetics of the spin spirals. 

5.1 Formalism 

SDFT involving c?-electron spins reduces approximately to a Stoner model l86l l87l . This reduction has 
the conceptual advantage of cutting the SDFT self-consistency loop into a band-structure part which for a 
given site and orbital-dependent exchange splitting, A, yields the site and orbital-dependent spin-moment, 
m (A) , plus a self-consistency condition which simply states that A = ml, where / is the Stoner (~ 
Hund's rule) interaction parameter. The band-structure part gives insight into the spin response of the 
non-interacting system, and not only in the linear regime. 

The spin arrangements which we shall consider are simple spin spirals. For these, the moment lies in 
the (x, 2/)-plane and has a constant magnitude, but rotates from site to site by an angle, cp (t) = q • t, 
proportional to the projection of the lattice translation, t in Eq. (Q~|i, onto the spin-spiral wave vector, q. 
Hence, the spin spiral with q at f produces FM order and the one with q at Y produces stripe order 
because the moment rotates by it upon y-translation, and by upon x-translation. Finally, the spin spiral 
with q at M produces checkerboard order because the moment rotates by tt upon y as well as upon x- 
translation. These spin spirals with q at high-symmetry points are collinear and commensurate but with 
the formalism which we now explain any q can be treated. 

In order to solve the band-structure problem in the presence of such a spin spiral, we use a basis set of 
localized Wannier orbitals times pure spin-functions, |j) and \l) , with quantization direction chosen along 
the local direction of the moment. In this representation, the one-electron Hamiltonian is translationally 
invariant, albeit with q-dependent hopping integrals, so that there is no coupling between Bloch sums with 
different wave vectors. As a consequence, the band-structure problem can be solved for any q without 
increasing the size of the primitive cell. [88 1 When merely seeking insights in this section, we shall neglect 
the interlayer coupling and use the 2D bands in the large BZ. So in this case, configuration space is in- 
variant to the t-mirror group, and spin space is invariant to the t-spinrotation group, which both have the 
same irreducible representations. As long as spin and orbital spaces remain uncoupled (spin-orbit coupling 
neglected), the one-electron Hamiltonian therefore factorizes down to the orbital and spin degrees of free- 
dom for a primitive cell of the t-group. This, together with SDFT, enables simple calculation of spin-spiral 
band structures, moments, and magnetic energies. 

The Hamiltonian turns out to be simply: 



(k) = 




in the local (y, I) representation and with the origin of k shifted to q/2. If the two paramagnetic Hamilto- 
nians, h (k) and h (k + q) , are identical (not merely their eigenvalues), this form is block diagonal. This 
is also the form appropriate for A larger than the bandwidths. For small A, it is more practical to transform 
to the (t T I) / V% representation in which 



In these expressions, h (k) is the paramagnetic 8 x 8 pd Hamiltonian whose eigenvalues, e a (k) , are the 2D 
bands discussed in the previous section, and A is an 8 x 8 diagonal matrix whose diagonal elements have 
the same value, A, for all five Fe d orbitals, and for all three As p orbitals. The approximation that only 
like orbitals couple goes back to the assumption that the spin density on Fe is spherically symmetric, and 
it is justified by the fact that, using this simple form for A, we find good agreement with the results of full 
calculations using SDFT for the spin-polarized bands. Of course, a form with the proper point symmetry 
on Fe could be used, but that would require more parameters. 

Although it takes the spin-spiral representation to see that the Hamiltonian above is general, we do 
note that, for an ^-electron in a commensurate antiferromagnet, a Hamiltonian of the form (TT2b is obtained 




(13) 
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by elementary means using the (ft, J|)-sublattice representation, global spin directions, and purely spatial 
Bloch functions. The Hamiltonian obtained for a ^-electron is the same, but with ft and JJ. interchanged. 

Diagonalization of the 16 x 16 Hamiltonian ( fT3T > yields energy bands, ep (k), and corresponding 
eigenvectors, {cft-i)i,p (k) , C(+ + n^p (k)} = {c^p (k) , ci t p (k + q)} , with I enumerating the 8 orbitals. 
(Here, ci : p (k) is a simplified notation for one of the 16 eigenvector components; for small A, this equals 
one of the 8 eigenvector components of the paramagnetic Hamiltonian, times 1 / \/2). We can now find the 
orbital-projected spin polarization of state (3k as: 

p^(k) = 2^(k)c^(k + q), (14) 

and summing this over the Fe d orbitals, pp (k) = J2i=i 5Pl,P s ) > an< ^ over tne occupied -or empty- 
states, we obtain the Fe moment: 

w — ,occ w — »cmpty 

-( A ) = E, k ^( k ) = -E /3k ^( k )- ( 15 ) 

This is the magnetic output of the spin-spiral band-structure calculation. 

5.2 Stripe 2D band structure 

In order to demonstrate how the spin-spiral formalism works for the 2D band structure of LaOFeAs, we 
start from the paramagnetic bands, e a (k) , decorated with the weight of each of the eight Wannier orbitals 
in Fig. [7j We consider a Y stripe, and thus prepare for the A-coupling as shown in the upper half of Fig. 
[T4l On top of the bands at k (in green) we place those at k + ny (in grey). Specifically, on top of the green 
fX bands we place the grey YM bands (which are the same as the XM bands with xz and yz exchanged), 
on top of the green XM bands we place the grey MX bands, and on top of the green f Y bands (which 
equal the fX bands with xz and yz exchanged) we place the grey YT band. The Mr bands couple with 
the XY bands, but since the latter were not shown in Fig. [7j this line is not shown in Fig. [T4l either. 

Now, the effect of A is to split degeneracies, e a (k) = ep (k + q) , by A times the geometrical average 
of the d characters, J2i=i 5 c i a s ) ci.p (k + q) . This of course only holds as long as A is so small that no 
further bands get involved. States without common <i-character therefore do not split. We note that states 
throughout the band structure split, independent of the position of the Fermi level, i.e. of the doping, but 
for small A only those states which gap around the Fermi level contribute to the magnetization and the 
magnetic energy, so this is how doping enters. 

The paramagnetic bands are seen to be linear inside an energy window of ±0.1 eV, at the most, around 
the Fermi level, and this means that effects of the exchange potential ±i A can be treated with linear- 
response theory only when A < 0.2 eV. 

5.2. 1 Bands and Fermi surface in the linear-response region, A < 0.2 eV 

A close inspection of the top xy panel of Fig. [I4]reveals that the crossing of the purely xy-like band with 
itself halfway between T and Y is 0.3 eV below the undoped Fermi level and thus requires A > 0.6 eV 
to gap around e p . On the other hand, the crossing of the xy hole band along FY with the yz jx electron 
band along Yf occurs only slightly below e p , meaning that the f -centered hole pocket and the Y-centered 
electron superellipse almost nest along f Y. This can be seen in Fig.[l5]where we show the Y-folded Fermi 
surfaces in brown lines. However, the yz/x band has only very weak xy character caused by its weak 
hybridization with the below-lying xy band, as was explained in connection with Eq. (f2]). Hence, due to 
lack of common orbital characters, these two states gap by much less that A. Finally, the xy hole band 
along fX crosses the xy/z electron band along YM at ~0.1 eV below ep, but due to the reduced xy 
character of the xy/z band, A must exceed ^0.3 eV to gap that part of the Fermi surface. As a result, for 
the value A=0.18 eV which via Eq. ( TTSb produces the same moment as the one observed experimentally, 
m (0.18 eV) = 0.3 \ib /Fe, the xy hole pocket does not gap. The Fermi surface calculated for A=0. 18 eV 
is shown by black lines in Fig. [15] For the T and Y-centered sheets, it only differs from the one calculated 
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Fig. 14 Top: 2D paramagnetic bands decorated like in Fig. [7] and prepared (folded) for stripe order with q=7ry : 
On top of the FY bands in green we place the YT bands in grey, on top of the TX bands in green we place the 
YM bands in grey, and on top of the XM bands in green we place the MX bands in grey. Bottom: 2D stripe band 
structure decorated with the orbital-projected spin-polarizations as given by Eq. ( 114b and with positive and negative 
polarizations in respectively dark and light blue. For the exchange potential, the value A=1.8 eV was used, which 
by Eq. J 1 5b yields the moment m(1.8eV) = 2.2 ps/Fe and corrresponds to the value 7=0.82 eV of the SDFT 
Stoner parameter. The dashed line is the Fermi level, which has moved up by 0.5 eV. Note that the paramagnetic 
and spin-spiral band structures are lined up with respect to the common paramagnetic potential, i.e. h (k) in the TB 
Stoner calculation. The 2D stripe Fermi surface is shown in Fig.[l6] For stripe order, As p projections cannot be 
spin-polarized and have therefore been omitted. 



for A=0 and shown in brown lines, because the Fermi level is slightly shifted due to gapping of the other 
sheets, i.e. those centered at M and X. That gapping, which causes the small moment of 0.3 /is/Fe, takes 
place as follows: 

The side of the X-centered electron superellipse which is normal to the x direction, i.e. which points 
towards T, matches the inner, longitudinal M-centered hole pocket both in Fermi-surface dimension (nest- 
ing) and in orbital character, xz. Those two bands therefore gap around the Fermi level, while the outer, 
transversal ( yz) hole sheet stays intact. Near the X and Y directions, the X electron sheet however matches 
the outer, transversal hole sheet in size and orbital character, Yz and Xz, respectively. So near those di- 
rections, the outer, transversal hole sheet is gapped while the inner, longitudinal sheet is intact. Finally, 
due to lack of common characters near the y direction, where the xy/z electron band does not hybridize 
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Fig. 15 Right: Nesting of orbital-projected Fermi-surface sheets for undoped 2D LaOFeAs (from Fig. [8} for stripe 
order with q=7ry (green arrow). Left: Partly gapped Fermi surface resulting from the small exchange potential 
A=0.18 eV which yields the small moment m=0.31 [lb /Fe. This corresponds to the Stoner parameter 1=0.59 eV. 

with the lower-lying xz band, no gapping occurs. As a consequence, small paramagnetic electron pockets 
with xyjz and transversal, xz/XY/zz/x characters occur. Such electron pockets will remain at the Fermi 
level, also for large A, as we shall see below. Note that the FS parts not gapped away are essentially not 
spin-polarized; this will not be the case for larger A. Since the Y stripe has antiferromagnetic order in the 
y direction and ferromagnetic order in the x direction one might expect higher conductivity in the x than 
in the y direction. However, most of the FS has a predominantly y-directed group velocity and predomi- 
nantly j/z-electron or yz longitudinal-hole character. Where the velocity is in the x direction, the character 
is xy/z electron or yz transverse hole. The xy hole pocket is isotropic in the plane. In conclusion, the 
exchange potential needed (A=0.18 eV) to give the observed moment with the Stoner model, is smaller 
than the fine structure of the bands. The gapping of the Fermi surface and the susceptibility, m (A) / A, 
therefore depend crucially on the details of the \t-and-orbital nesting. 

5.2.2 Bands and Fermi surface beyond the linear-response region, 0.2 cV < A < 3 eV 

The exchange potential obtained selfconsistently from the SDFT and yielding the proper crystal structure, 
is ten times larger: A=1.8 eV, and thereby has the same scale as the structure of the subbands, i.e. this A 
is intermediate and linear-response theory invalid. The SDFT value of the Stoner parameter is 7=0.82 eV 
and the moment is m (1.8 eV) = 2.2 /Us/Fe. Compared with the maximum moment of 4 ns/Fe for the d e 
configuration, the value 2.2 /is /Fc is intermediate. The lower part of Fig. [l4]now shows the stripe bands 
for this situation. These bands are complicated, because the gapping and spin-polarization depend on the 
energies and the p and d orbital characters of those bands at k and k + q which are separated by less than 
~ A. In the present section, we shall describe those bands and their Fermi surface and calculate specific, 
important levels analytically. 

In Fig. Q31 the paramagnetic and the spin-spiral band structures are lined up with respect to the common 
paramagnetic potential, i.e. h (k) in the TB Stoner calculation. We see, as was pointed out before, that 
bands with d character split irrespective of their position relatively to the Fermi level, but those in the 
lower half the d-band structure generally shift downwards (dark blue) while those in the upper generally 
half shift upwards (light blue). We recall that the shift upon an increase of the exchange potential is the 
negative of the spin-polarization: dep (k)/ d (A/2) = — p@ (k) , by lst-order perturbation theory. In fact, 
there happens to be a fairly well-defined dividing line between positively and negatively spin-polarized 
bands around eV. Moreover, on this dividing line, the non-hybridizing yz and zz/XY bands along 
FY are nearly degenerate and dispersionless, so increasing A beyond 1.8 eV will open up a gap which 
extends throughout the BZ and makes any correponding d 5 material (e.g. LaOMnAs or LaOFeN (89)) an 
antiferromagnetic insulator. 

Copyright line will be provided by the publisher 



32 



O. K. Andersen: On the multi-orbital band structure and itinerant magnetism of iron-based superconductors 




Fig. 16 Fermi surface for stripe order (q=7ry) resulting from the large exchange potential A=1.8 eV which yields 
the large moment m—2.2 fis /Fe corresponding to to the SDFT Stoner value J=0.82 eV. See also Fig,sll5land[8l 



We now discuss the intermediate-moment stripe bands for 2D LaOFeAs in detail. Starting again on the 
left-hand side of Fig. [14] with the crossing of the paramagnetic, nearly pure xy band with itself, halfway 
between T and Y at —0.3 eV, we see the bands split to the energies —0.3 ± 0.9 cV = — 0.3eV ± A/2 
around the Fermi level, which has now moved up to 0.5 eV, with the lower and upper bands fully spin- 
polarized. This gap extends in a large region around the TY-line. So, whereas for small A, the T xy hole 
and Y xy-yz electron sheets did not gap at all, for intermediate A, these two FS sheets no longer exist. It 
is however only the xy-parts which are gapped away: Due to strong yz/x hybridization, the yz-like bands 
near f (green) has no partner at Y (grey) within the ±A/2-range with which it can couple. This band 
therefore only splits midways between f and Y, i.e. near ^Y, but hardly closer to f and towards X (grey). 
Hence, the reason for the disappearance of the yz-pait of the superellipse at Y is not gapping, but the 0.5 
eV upwards shift of the Fermi level. This shift is due to the fact that with configuration d 6 , the Fermi level 
lies in the upper half of the rf-like band where most bands are shifted upwards by the exchange potential 
and drag the Fermi level along with them. The shift deF (A)/ d (A/2) is upwards, if at ef (A) the density 
of 4- states exceeds that of j" states. It may be noted that the crossing of the grey xy/z and yz bands along 
YM, which is at d 6 for the paramagnetic bands, still occurs for the stripe bands, but far below the Fermi 
energy. 

The bands expected to gap mostly for Y-stripe order are d bands dispersing less than A/2 in the y 
direction, i.e. the xz and zz bands. The xz band is dispersionless near the XM line where it forms the 
M-centered transverse hole band and the bottom of the X-centered electron band. Further towards the TY 
line, however, the fc a -dispersion of the xz band becomes large due to by-mixing of y character towards T, 
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and also the concomitant dilution of d character reduces the exchange coupling. As a consequence, near 
the XM line, the xz electron and hole bands split to ~ ±A/2, whereby the 4- band is above the shifted 
Fermi level and the j" band is far below. This emptying of antibonding xz states yields the observed 0.5 
% orthorhomic contraction in the x direction, along which the spins are aligned ferromagnetically ll37ll . In 
the remainder of the zone, there are no bands at the Fermi level. The paramagnetic zz band has a width 
of about 1 cV ~ A/2 and is centered at ~ —0.6 eV, so we expect the exchange splitting to shift the zz I 
band up to -or above- the shifted Fermi level. Fig. [14] shows that this is roughly the case, but the details 
are more complicated, due to strong hybridization with the XY band, as we shall see later in connection 
with Eqs.dPTli and d 1 8t . The paramagnetic and stripe bands, albeit merely for positive energies and the 
exchange splitting, A=l.l eV, may be seen more clearly in the left-hand side of Fig. [18] We shall return to 
this figure. 

Only two bands remain at the Fermi level: (1) the longitudinal yz I band dispersing downwards from 
M, crossed by and hybridizing with a weakly spin-polarized xy/z j" band dispersing upwards from X, and 
(2) the zz l band dispersing downwards from M, hybridized with more weakly spin-polarized XY j" band. 
Being respectively even and odd by reflection in a vertical mirror containing nearest-neighbor As atoms 
and the q vector (see Fig. [2}, bands 1 and 2 cannot hybridize and therefore cross at a Dirac point along 
the XM line. This pins the <i 6 -Fermi level and gives rise to a Fermi surface shaped like a propeller, |[T8l 
with two electron blades and a hole hub (Fig. IT6b . The hub is yz I like and the inner parts of the blades 
are mixed zz I - XY f like, while the outer edges are mixed yz 4- - xy/z f like. Due to the pinning, the 
propeller shape is robust, e.g. not sensitive to A. Hole doping by a few per cent will bring the Fermi level 
to the Dirac point and make each electron blade shrink to a point. Upon further hole doping, the blades 
will reappear as hole sheets. The velocity of the xy/z j" and the zz I parts of the anisotropic cone are 
respectively ~ 1 eV • a = 2.9 eV A and ~ 0.4 cV • a = 1.1 eV A. These Dirac cones in the stripe-ordered 
SDW state have been predicted [78 1 and later observed using respectively quantum oscillations |3~T1 and 
ARPES f79l . Compared with ours, the experimental velocities seem to be renormalized by a factor ~ 1/4. 
This Dirac cone will be gapped by any lattice imperfection breaking the above-mentioned mirror symmetry 
and is therefore not "protected." 

In order to explain how the complicated spin and orbital characters of conduction bands 1 and 2 arise 
from the paramagnetic band structure, let us consider the simple case that the paramagnetic TB Hamilto- 
nian h (k) in Eq. ( fT~3T > is a 2 x 2 matrix. Its eigenvalues are the bonding, eb (k) , and antibonding, e a (k) , 
paramagnetic bands with the respective eigenvectors, {— sin (k) , cos</> (k)} and {cos<fi (k) , sin0(k)} = 
{c (k) , s (k)} (we have chosen the phases of the orbitals such that the Hamiltonian is real). Taking the 
first orbital is a d and the second as a p orbital, transformation of the spin-spiral Hamiltonian (fT~3l > to the 
bonding-antibonding representation yields: 



pd,q 



(k) = 



/ £a (k) 

e b (k) 

c.c. c.c. 

c.c. c.c. 



iAc(k)c(k + q) 



■|A«(k)c(k- 
e a (k + q) 




-iAc(k) S (k + q) \ 
|A s (k) s (k + q) 


Sb (k + q) / 



(16) 



because the pp and pd elements of the exchange block, ^A, vanish. Here, c 2 and s 2 = 1 — c? are the 
d characters of the antibonding and bonding levels, respectively. Note that only the d -but not the p— 
characters need to be the same at k and k + q. This 4x4 form ( TTol i is exact when k, and thereby k + 7ry, 
is at a high-symmetry point, because here, h (k) factorizes into blocks of dimension < 2 x 2, with the 
XY / zz/ x-b\ock. at X as the notable exception (see caption to Fig.[2]i. We thus start explaining the spin 
and orbital characters of band 1 by using ( fl~6b to couple the levels at k = X to those at k + q = M : 
the paramagnetic xy-like antibonding and bonding levels at X (see Figs.P71and[T4Tl are the strongly xy-\\ke 
bottom of the electron band (green), for which e a (X)=— 0.46 eV and c (X)=0.95, and the strongly y-like 
level at eb (X)=— 2.79eV. These levels couple to the xy-like levels at M (grey) of which the antibonding 
one at e a (M)=0.95 eV is mostly z-like, and the bonding one at Sb (MJ=— 3.36 eV is mostly xy-like, 
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s (M)=0.83. The four rry-like stripe levels at X are thus the eigenvalues of the Hamiltonian: 

( -0.46 (A/2) (0.95) (0.56) - (A/2) (0.95) (0.83) ^ 

-2.79 -(A/2) (0.31) (0.56) (A/2) (0.31) (0.83) 

c.c. c.c. 0.95 

V c.c. c.c. —3.36 J 

These eigenvalues (seen in Fig.fPfli are: 1.12, —0.45, —2.73, and — 3.59cV when A/2=0.9eV, and hence 
little perturbed by the stripe order. The reasons are that the paramagnetic levels are separated by more than 
A/2, and that the p hybridization reduces the geometrical averages of the d characters far beyond unity, 
except for two levels which are, however, separated by as much as 3 eV. In particular the state of interest, 
the one at —0.45 eV, has been pushed down by the level at 0.95 and up by the one at —3.36, both at M, 
and as a result, has moved by merely 0.01 eV. For the same reason, its spin polarization is only about 50%. 
This state thus remains essentially the (green) bottom of the electron band at X. 

The paramagnetic y z-like antibonding and bonding levels at X are the strongly y z-like bottom of the 
longitudinal hole band, for which e a (X)=— 1.26eV and c(X)=0.998, and the strongly z-like level at 
Eb (X)=— 2. 12 eV. These levels couple to the yz-like levels at M (grey), which are top of the doubly degen- 
erate hole band, for which e a (M)=0.21 eV and c (M) =0.90, and the y-like level at e b (M) =- 1 .97 e V. With 
these values, the y z-like eigenvalues of the 4x4 stripe Hamiltonian ( fToT l becomes 0.58, —1.39, —2.11, 
and —2.21 eV. Here the uppermost level, being near the Fermi level, is our band 1. It is described to a 
good approximation by using merely the antibonding paramagnetic states, i.e. by the 2x2 Hamiltonian 

s a (k) A c ( k )c(k + q) \ / -1.26 (A/2) (0.998) (0.90) \ 

|c(k)c(k+q) e Q (k + q) J \ c.c. 0.21 J 6V ' 

because these states are the only ones with substantial yz character, as may also be seen from Fig. [14] The 
uppermost state thus has about 80% Myz and 20% Xyz character, and its spin polarization is —75%. 

As we now move from X towards M, via |XM as in Figs. [14] and QjO the Xxy/z j" band disperses up- 
wards from —0.45 cV like the paramagnetic (green) electron band, and the yz I band disperses downwards 
from 0.58 eV like the paramagnetic (grey) longitudinal Myz hole band. At about 1/4 the distance to M, 
these two bands suffer an avoided crossing. From there on, the xy/z t band continues towards M like the 
paramagnetic band, apart from the facts (a) that it gets folded and split at |XM with a covalency-reduced 
A, and (b) that it continues to hybridize with the longitudinal yz l band whose downwards dispersion 
(grey) is halted near |XM due to repulsion from the upcoming (green) yz band. The xy/yz hybridization 
matrix element along XM, — 2\/2 (t^y Xz + t\lvz + ^xy Xz) s ^ n ^2/> would have vanished, if the xy and 
Xz Wannier orbitals of the pd set shown in Fig.|6]had been respectively symmetric and antisymmetric with 
respect to the Fe plane. Band 1 obtained from the pd Stoner model is somewhat more shallow than the one 
obtained from a standard LAPW calculation, which yields a bandwidth of 0.6 eV. This can be traced back 
to the gap between the downwards-dispersing longitudinal yz \, band and the upwards dispersing xy/z j" 
band, being 0.5 eV in the model but merely 0.2 eV in the LAPW calculation, thus causing the model band 
to be 0.3 eV more narrow. Reducing Xz + t** Xz + t^J j, does not entirely remove this discrepancy, 
which might also be due to our assumption of a spherically symmetric exchange potential. 

We now come to band 2. It was pointed out in Sect. 13. 31 and can clearly be seen in Figs. l7land[T4l 
that the paramagnetic XY and zz bands hybridize strongly, except along the f M lines, and have avoided 
crossings around X and Y causing them to gap around the Fermi level. They also gap around the Fermi 
level along TM, but due to avoided crossings with other bands. In order to understand the effect of a stripe 
potential, let us use a 4 x 4 model like ( fTST i, but now for two d orbitals. In this case, the exchange block is 
constant and the spin-spiral Hamiltonian in the bonding-antibonding representation becomes: 



Hdd^ (k) 



/e.(k) ±Acos(0(k + q)-0(k)) -±Asm(</)(k + q)-</>(k)) \ 

e 6 (k) iAsin(<£(k + q)-«Mk)) |Acos (0 (k + q) - (k)) 

c.c. c.c. e a (k + q) 

V c.c. c.c. Eb (k + q) I 
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(17) 

Clearly, if the bonding and antibonding linear combinations of the two d orbitals were the same at k and 
k + q, i.e. if <fi (k) =0 (k + q) , then the bonding-antibonding representation (T7\ would be identical 
with the dd representation. In that representation, the off-diagonal block is times the unit matrix 
because we have assumed the exchange potential to be spherical. If now also e a (k) =e a (k + q) and 
Eb (k) =Eb (k + q) , as would be the case at the zone boundary, then we could transform to the local spin 
representation (fT2] i and would then immediately realize that the 4 stripe eigenvalues (in eV) are: 

e a QxMjiy «0.8±0.9= j j and e b QxMJ ±y « -0.8±0.9 = 

(18) 

This fits well with Fig.[l4] from where we have taken the values 0.8 and —0.8 eV, for the paramag- 
netic antibonding and bonding levels at |XM. The stripe level belonging to band 2 is the bonding, 
minority-spin level at 0.1 eV. To be honest, the energies ±0.8 eV of the paramagnetic levels do include 
weak hybridizations with the xz and x orbitals, which have been neglected in the 4x4 (XY, zz) model 
( fTTI ). We should also warn that although h (k) and h (k + q) have the same eigenvalues at the zone 
boundary, these matrices are generally not identical; off-diagonal elements may have different signs, e.g. 
hxY.zz (-§ y) = h XY ,zz (f y) but h xx>zz (-§ y) = -h xz , zz (f y) . 

As seen from the figures, the paramagnetic bonding and antibonding XY/ zz bands disperse less than 
their separation, 5 ~ 1.6 eV, so in order to be able to diagonalize the 4x4 stripe Hamitonian (fTTI i let us 
stay with the assumption that e a (k) = e a (k + q) = 6/2 and Eb (k) = Eb (k + q) = —5/2, but drop the 
assumption that <p (k) = <\> (k + q) . The 4 stripe bands are then given by: 



e(k) = ±^y/5 2 ±2SAcosip(k) + A 2 , where tp (k) = (k + q) - <f> (k) . (19) 

For tp=0, the zone-boundary case d 1 3b . the bonding and antibonding d orbitals are identical and each of 
the four levels, ± (6/2) ± (A/2) , have pure spin and pure bonding or antibonding character. For </j=7r/2, 
the bonding and antibonding bands have orthogonal d characters and therefore only have off-diagonal 
exchange coupling. This does not split the spin-degenerate bonding and antibonding levels, but separates 
them to ±i yj 5 2 + A 2 . The realistic case at X (green) is that the antibonding and bonding levels have 
roughly the same XY and zz characters. We therefore take (j> (X)=7r/4, and the XY orbital as the first 
orbital. At M (grey) the "antibonding" level has pure XY and the "bonding" level pure zz character, so 
<p (M)=0. As a result, the bands have dispersed from the levels given by CUi to ±\\/5 2 ±x/2<5A + A 2 = 
±1.57 and ±0.66 eV at X=M. This agrees well with what is seen in Fig.[T4l and explains why band 2 
disperses upwards from 0.1 at |XM to 0.66 eV eV at X. The latter level is essentially bonding zz \. 
/XY t- 

Armed with the detailed understanding of the interlayer hopping provided in Sect.|4] and that of the 
generic 2D stripe band structure provided above, the interested reader should be able to digest the compli- 
cated 3D stripe bands for specific materials found in the literature and recently reviewed in Refs. (HID. 

5.3 Magnetization and magnetic energy 

Apart from the spin-spiral band structure described above in the case of stripe order, the output of a band- 
structure calculation with an imposed exchange potential, A, is the Fe-magnetization, m (A). In Fig. [17] 
we now give the results obtained for stripe (q = 7ry) and checkerboard (q = 7rx + 7ry) orders for various 
electron dopings, x, in the rigid-band approximation. With the Stoner approximation, we have been able 
to afford sampling the spin-polarization over a very fine fc-mesh so that nesting features are resolved. 
Since the magnetization increases linearly with A, when it is small, we plot the static spin suceptibility, 
X (m) = m (A) /A, and as a function of m rather than of A, since the m is an observable. Note that with 
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Fig. 17 Non-interacting, static spin susceptibilities \ (m) = m/A for 2D LaOFeAs calculated from spin-spiral 
band- structure calculations, i.e. from diagonalizing Hq (k) in Eq. d 1 3t and finding m from the eigenvectors according 
to Eq. dl5t . The paramagnetic TB pd Hamiltonian, h (k) , was calculated for the observed structure. The electron 
dopings, x (in e/Fe), were varied in the rigid-band approximation. For a given value of the Stoner interaction parameter, 
/, the self-consistent moment is given by \ (m) = 1/1 . The SDFT value of I is 0.82 eV and the value fitting the 
experimental moment and its doping dependence is 0.59 eV. 



4 — x empty bands, 4 — x is the value of the saturation magnetization and \ (m) therefore vanishes for m 
larger than this. 

Given a value of the Stoner exchange-coupling constant, /, the self-consistent value of the magnetization 
is the solution of the equation \ (m) = 1 //, and we see that for the SDFT value, 7=0.82 eV, m~2.2 /zs/Fe 
for both stripe and checkerboard order, and that m decreases with electron doping. The reason for the latter 
can be understood by considering the stripe bands for A= 1 .8 eV at the bottom of Fig. [14] The magnetization 
along the upwards-sloping line \ = m / (1-8 cV) in Fig. [I7]is the sum over the empty bands of their spin 
polarizations, taken with the opposite sign according to Eq (15[ , i.e. of the light-blue fatness. Since light- 
blue is seen to dominate over dark-blue, at every energy above the Fermi level, moving the latter up, as 
electron doping will do in the rigid-band spproximation, must decrease the moment. 

Coming now to the small-moment part of Fig. [T7] the linear response, \ (0) i i s seen to be particularly 
large for stripe order and no doping. This is due to the good nesting shown in Fig.[l5]between the xz part 
of the X-centered electron superellipse and those of the M-centered hole pockets. This nesting is, however, 
sensitive to the relative sizes of electron and hole sheets and is therefore rapidly destroyed with electron (or 
hole) doping, thus causing \ (0) to decrease. Also, increasing A beyond 0.2/1.8 ~ 0.1 eV is seen to make 
X (m) decrease rapidly. We should remember (Fig. H} that the top of the f -centered xy-like hole pocket is 
merely 0.06 eV above the Fermi level for the pure material and that this pocket disappears once the doping 
exceeds 0.1 e/Fe. We recall also, that the top of the M-centered hole pockets is merely 0.2 eV above 
the pure Fermi level and that these pockets disappear, as well, once the electron doping exceeds 0.3 e/Fe. 
The T-centered hole pocket and the .xy-part of the Y-centered electron superellipse start to gap when A 
exceeds 0.2 eV, and for larger A, the xy moment becomes as large as the xz moment. The experimentally 
observed moment in LaOFeAs is ~0.4^s/Fe, stripe ordered, and vanishes for x > 3%. This would be 
consistent with our 2D bands and the Stoner model if 7=0.59 cV. However, only by virtue of its large 
moment, ^2/i^, does the SDFT yield the observed large value of the As height, rpO.93, and the observed 
0.5% orthorhombic contraction in the direction of ferromagnetic order. 

For checkerboard order, q is at the M point and this places the T and M-centered hole sheets -as well 
as the X and Y-centered electron sheets- on top of each other. Nesting of sheets with the same (electron 
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Fig. 18 Band-resolved magnetic energy for stripe order (q at Y). Red: Unoccupied part of the folded paramagnetic 
band structure, i.e. e a (k) and e a (k + q) for A=0. This is the upper part of Fig. [14] but only for positive energies and 
without fatness. Black left: Band structure, ep (A,k) , obtained by applying an exchange potential of intermediate 
strength, A=l.l eV, yielding a magnetization of 1.5 /is /Fe and thus corresponding to 7=0.73 eV. The Fermi level 
(dashed line), ef (A) , has moved up by 0.3 eV. The dominant d-orbital characters have been written onto the bands. 
Black right: Unoccupied part of the band structure, ep (A, k) = ep (A, k) + jpp (A, k) A, corrected for double 
counting such that the magnetic energy gain per /?k-hole is ep (A, k) — ep (0, k) , i.e. black minus red. Note that only 
the unoccupied part of the bands are shown and, in particular, that the double-counting corrected e-bands are truncated 
by the e-band Fermi level, ef (A) , a truncation which does not occur at the same energy for all e-bands. 



or hole) character is not optimal for gapping, and \ (0) is therefore neither very high nor very doping 
dependent. 

For moments so low that for all possible spin orientations m(A) is linear and the magnetic energy 
quadratic, the effective coupling between the spins can be expressed in terms of the (Stoner enhanced) 
linear, static spin susceptibility. For moments so large that the system is insulating, on the other hand, 
the electronic degrees of freedom can be integrated out, whereby the coupling between the spins is given 
by a Heisenberg model. That model, with 1st and 2nd-nearest neighbor antiferro-magnetic couplings and 
J\ 2 J2 has, in fact, often been used to describe the magnetism of the iron pnictides. This is, however, 
hardly justified because the iron pnictides are metals, presumably with intermediate moments. Full SDFT 
calculations, like the spin-spiral calculation described at the end of this section, do however account well 
for many experimental observations and as a first step towards deriving better exchange models, we shall 
therefore try to explain the origin of the magnetic energies using the Stoner model. 

The relation between the magnetic energy for a particular spin spiral, a Y stripe, and the underlying band 
structure is illustrated in Fig. [18] Its left-hand side shows the paramagnetic bands for positive energies, 
i.e. the unoccupied bands, as well as the stripe bands for the somewhat reduced value A=l.l eV of the 
exchange splitting. This corresponds to 7=0.73 eV and to an SDFT calculation for BaFe2As2 adjusted 
to the experimental dHvA FS ||6T1 . The dashed line shows the Fermi level, ep (A)=0.3 eV. Note that, 
like in Fig.[l4l the paramagnetic and spin-spiral band structures are lined up with respect to the common 
paramagnetic potential. We clearly see that the A=l.l eV stripe perturbs all bands, but the p-like ones the 
least. 

In SDFT, the total energy is a stationary functional of the electron and spin densities. For densities 
which can be generated by occupying the solutions of a single-particle Schrodinger equation for a local 
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potential according to Fermi-Dirac statistics, the value of this functional is simply the sum of the occupied 
single-particle energies, minus corrections for double-counting of the Hartree and exchange-correlation 
energies. This holds when the potential is the self-consistent one, i.e. the one which minimizes the energy 
functional. For our Stoner model with h (k) describing the self-consistent paramagnetic bands and for A 
taking the self-consistent value, m (A) I, the double-counting correction of the magnetic energy is simply 
jm (A) A = Y^ik \p@ ( A, k) A. For d 6 materials, we shall sum over empty states, because of those there 
are only 4/Fe, and nearly all have negative spin polarization. The double-counting correction of the stripe 
bands has now been performed on the right-hand side of Fig.[T8]from where we realize that these (black) 
bands, eg (A, k) = eg (A, k) + jpp (A, k) A, are/ar less perturbed than the real stripe bands, ep (A, k) . 
This means, that the state-resolved magnetic energy gain (black minus red), ep (A, k) — Ep (0, k) , is 
concentrated near the exchange gaps and near the paramagnetic and spin-spiral Fermi surfaces. Note that 
each of these empty bands, ep (A, k) and Ep (0, k) , should be defined as = Ep (0) , if that band is 
occupied. This means that the empty paramagnetic bands, ep (0, k) , are continuous, but truncated with 
a kink at the lower figure frame. The empty, corrected magnetic bands, ep (A, k) , should be truncated 
discontinuously. Due to the 0.3 eV upwards shift of ef (A) , all empty parts of the corrected bands are 
above the frame of the figure so that all empty magnetic bands are visible. 

In order to see which states contribute to the energy of the Y stripe, let us now once again start from |Y 
and move to the right in the right-hand figure. The zone-boundary gapping around 1.2 eV of the corrected 
XY/zz/xz/y band is small and fairly localized near iY, and there, its positive and negative contributions 
nearly cancel. The gapping around 0.7 eV of the yz/xy/x band is much larger due to the dominating yz 
character of this band, but here again, the negative and positive contributions essentially cancel, until k 
gets closer to f . There, the character of the lower band is xy from the f -centered hole band and xy/z from 
the Y-centered superellipse electron band. This band contributes positively to the energy of the stripe in 
the large region of the BZ around f =Y where the FS is completely gapped, a contribution which integrates 
up to about half the stripe energy. Between T and X, the black xy/z band is seen to suffer an avoided 
crossing with the band formed from the X-centered xz jy electron band and the longitudinal M-centered 
xz hole band. Near the avoided crossing the magnetic energy density becomes negative, but is essentially 
cancelled by the positive contribution from the upper band. Closer to X, the contribution from the lower 
band becomes positive again, and its positive magnetic energy density is seen to extend over the large 
region around X=M where the xz part of the FS is completely gapped. This part of the band is formed 
by coupling of the paramagnetic, flat X-centered xz electron band, which is occupied and therefore lies 
below the frame af the figure, and the paramagnetic transversal M-centered xz hole band, which becomes 
occupied outside the transversal hole pocket. So whereas the empty, black xy-xz band extends smoothly 
throughout the BZ, its paramagnetic partner, against which we measure the band-resolved magnetic energy, 
goes to zero at a few places in the BZ, such as outside the transversal M-centered hole band along XM. 
The magnetic one-electron energy of the xz-\ike band thus reaches +0.5 eV between X and |XM. 

Also the two uppermost black bands, which are degenerate at T and then split into longitudinal and 
transversal p/t bands, stay empty, i.e. they extend smoothly throughout the BZ. The lower of these bands, 
the one which is longitudinal p/i-like near f , becomes M z/xy like near X and is seen to contribute 
negligibly to the magnetic energy in LaOFeAs (but possibly not in the bet materials). The same holds 
between |XM and ~ -|XM for its lower partner, which is essentially the paramagnetic X-centered xy/z 
electron band extending in the direction towards M. At |XM, where the magnetic band becomes the tip of 
the propeller blade, the corrected xy/z band jumbs to 0, whereby the one-electron magnetic energy jumps 
from to — 0.25eV. Also the M-centered yz hole band coincides with the corrected magnetic yz band, so 
that only the Fermi-surface truncations contribute to the magnetic energy, which in this case amounts to a 
small negative energy from the region between the propeller hub and the yz-part of the M-centered hole 
pockets. 

The higher of those two bands which are degenerate at T, i.e. the one which is transversal p/t near T, 
develops into the uppermost of the four bands formed by the coupling of the paramagnetic X XY/ zz and 
MXY bands described in connection with Eq. ( fT9l l, albeit for a larger A. This band is seen to contribute 
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Fig. 19 k-resolved magnetic Y-stripe energy, ~^2 c ™ pty [ep (A, k) — e$ (0, k)] , as in Fig.Q~8]but for A= 2.3 eV and 
summed over the empty bands and shown throughout (k x , fcj / )-space. Magnetic energy gains are red and losses blue. 

The discontinuities caused by the propeller sheet are clearly seen. While the left-hand figure is for 2D LaOFeAs 

i i 

with the observed structure (ry=0.93) , the right-hand figure shows the result obtained by reducing ti£ z from 0.52 to 
0.30 eV. With the latter value, the 2D paramagnetic z/xy band becomes degenerate with the top of the hole bands at 
M, and thus forms a Dirac cone with the longitudinal t/p band. This corresponds to the elongation rj^l.2. 



positively to the magnetic energy, a contribution which is, however, overwhelmed by a large, negative 
contribution from the second of the four bands, which is part of the XY/zz band decreasing from 1.3 eV 
at ^Y to 0.6 eV at X, and finally to 0.2 eV around |;XM. At this point the zz/XY band is on the inner 
part of the propeller blade, and therefore jumps discontinuously to 0, where is stays until reaching the zone 
boundary at |XM. The magnetic one-electron energy of the zz/ XY-ltite band thus jumps from —0.6 to 
-0.8 eV at |XM. 

It is simpler to eyeball the balance of magnetic one-electron energies along the line between X and 
|XM, where we found large cancellations, if we connect the magnetic XY/zz and xz bands according 
to energy. In fact, the real zz/XY and xz bands may cross along f X and YM, but not between X and 
iXM. When doing so, we see that the magnetic energy loss from the XY/zz-xz band lying near 0.5 eV 
is nearly balanced by the gain provided by the upper XY/zz band lying near 1 eV. The lower xz-zz/XY 
band lies ^0.2 eV above the xz part of the paramagnetic M-centered hole band and thus gains considerable 
magnetic energy. The Fermi-surface contribution seems to be small because kp y for the outer, transversal 
M-centered hole surface is about the same dimension as the distance from the center of the hub to the inner, 
zz/XY part of the blade. Finally, by joining this xz-zz/XY band across the blade to the xy/z band, we 
see that there is a loss of one-electron magnetic energy of about 0.2 eV inside the blade. 

This is seen quite clearly in the left-hand part of Fig.Qjlwhere we show the k-resolved magnetic energy, 
which is the state-resolved magnetic energy considered above, summed over empty bands. What stabilizes 
stripe order when its moment is > 1 ^b/Fc, is then, first of all, coupling of the paramagnetic F-centered 
xy hole and the Y-centered xy/z electron bands over a large part of k-space centered at f . The second, 
almost as large contribution comes from coupling of the xz part of the paramagnetic M-centered hole band 
to that of the X-centered electron band over a smaller part of the zone centered at X. These T and X- 
centered red regions do not overlap. The Fermi-surface contributions to the magnetic energy are relatively 
small, and the positive (red) contribution from the hub tends to cancel the negative (blue) contribution from 
the blades. This being the case, it should be possible to derive a Heisenberg model which fairly accurately 
describes the change of the magnetic energy for pertubations of q around commensurable stripe order. 

We can now address the interplay between the distance between the As and Fe sheets and the stripe 
magnetism: The main effect of increasing this distance, rj, is to decrease the z-xy hybridization, as was 
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Fig. 20 Magnetic moments (upper panels) and energies (lower panels) per Fe of spin spirals as functions of q for 
different electron (x) and hole (y) dopings in the virtual-crystal approximation. These results were obtained by self- 
consistent SDFT-LMTO calculations. |48J The energies of the Ji, J2 Heisenberg model for x(y) = and 0.2 are given 
by respectively dashed and dash-dotted lines. Representative real-space spin structures are shown at the bottom right 
for the q-vectors denoted by dots (adapted from Ref. l48l ). 



discussed in Sect. I4.ll and thereby to decrease the splitting between the paramagnetic z/xy antibonding 
and xy/z bonding levels near M, and thus to move the empty M z/xy level down. The z/xy electron 
band from Y=f to M=X, as well as the ones folded from X to |XM (see FigJT4t. will be less steep 
and more xy-like upon increasing This, in turn, will increase the polarization and the gapping of the 
xy-like stripe bands and thereby increase the fixed-A moment, m (A) , and the differential suceptibility, 
X = dm (A) / dA. The increase of the self-consistent moment will finally be enhanced over the increase of 
m (A) by the Stoner factor (1 — I\) ■ At the same time as the moment increases, so does the gapping, 
the stripe energy, and, hence, the magnetic energy. Also the mere flattening of the z /xy bands increases 
the magnetic energy by extending the regions around T and X of positive magnetic one-electron energy. 
As seen in the right-hand part of Fig. [19] the f region increases in the k x direction and the X region in the 
k y direction. This explains why spin polarization tends to increase the vertial Fe-As distance. 

In conclusion, what stabilizes stripe order when its moment is > 1 hb /Fe, is -first of all- the coupling 
of the paramagnetic f -centered xy hole and Y-centered xy/z electron bands over a large part of k-space 
centered at f . The second, almost as large contribution comes from coupling of the xz part of the para- 
magnetic M-centered hole band to that of the X-centered electron band over a smaller part of the zone 
centered at X. These T and X-centered regions do not overlap. The Fermi-surface contributions to the 
magnetic energy are comparatively small and mostly negative. This being the case, it should be possible 
to derive fairly accurate Heisenberg models describing the change of the magnetic energy for pertubations 
of q around commensurable stripe order. For a start, one might compute the spin-spiral energy dispersions 
as a function of q, using the simple TB Stoner model ( TT3T > and then analyse which one-electron states are 
reponsible for the energy changes, like we did for q=7ry. 

The results shown so far were obtained using the spherical Stoner model, which allowed us to simplify 
the calculation of the spin-spiral band structures and magnetic energies so much, that we might understand 
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the results by solving simple analytical problems. Although approximate, this model is in many respects 
more general than SDFT calculations. Now, coming to the end of our tutorial paper, we show in Fig. 
|20]results of SDFT spin-spiral calculations for LaOi_ x F x FeAs and Bai_2 y K2 y Fe2As2 of self-consistent 
moments and energies as funtions of q and doping in the virtual-crystal approximation (VCA). 11481 Note 
that in this figure, q takes the usual rXMT path, so that the spin-spiral patterns near X correspond to an X 
stripe. For LaOFeAs, the moment is seen to be ~^1.3 ^s/Fe for both stripe and checkerboard order, which 
is somewhat smaller than what we obtained in Fig.[l7j presumably because the moment is calculated by 
integration of the spin-polarization in an Fe sphere with radius ~a/2 in the SDFT calculation, rather than 
being summed over Fe Wannier orbitals. Other causes could be our use of the Stoner aproximation with 
too high an / and a spherical A. More significant is, however, that whereas we found, and understood, 
that in the rigid-band approximation the large moment decreases with electron doping, and increases with 
hole doping, the behaviour seems to be the opposite in the VCA approximation, both for electron-doped 
LaOi_ x F x FeAs and hole-doped Bai_2 y K2 y Fe2As2. The VCA approximates Oi_ x F x (or Bai_2 y K2 y ) by 
a virtual atom having a non-integer number of protons. Such anomalies have recently been discussed for Co 
and Ni-substitution at the Fe site using supercell calculations,] 90 1 but not for substitution in the blocking 
layers. We may speculate that at least for Bai_2 y K2 y Fe2As2, the strong Ba 5c? hybrididization near the 
Fermi level found in Sect. l4.2l could make substitution of Ba by K a non-trivial doping. Nevertheless, the 
spin-spiral energies shown in the lower half of Fig.[20]agree with experiments to the extent that the stable 
spin order is the stripe for low doping, but shifts to in-commensurable order for higher doping, x> 5% 
in LaOFeAs. This is consistent with the onset of superconductivity in this material, but is a completely 
different, and presumably more accurate scenario than the one suggested by Fig. [17] The black dashed 
and dotted lines are fits by a simple Heisenberg nearest and next-nearest neighbor (Ji, J2) model to the 
calculated spin-spiral dispersions for respectively the undoped and 20% doped compounds. These fits are 
not bad, although exchange interactions of longer range are needed to fit the LaOFeAs incommensurability. 
However, the SDFT energies of other spin arrangements, such as starting from stripe order and then rotating 
the spins on one Fe sublattice rigidly with respect to those on the other, could not be reproduced by the 
(Jl, J2) model. 



Summary 

We hope to have given a pedagogical, self-contained description of the seemingly complicated multi- 
orbital band structures of the new iron pnictide and chalcogenide superconductors. First, we derived a 
generic Fe d As p TB Hamiltonian by NMTO downfolding of the DFT band structure of LaOFeAs for the 
observed crystal structure. By use of the glide mirror symmetry of a single FeAs layer, its primitive cell 
was reduced to one FeAs unit, i.e. the TB pd Hamiltonian, h (k) , is an 8 x 8 matrix whose converged 
analytical expressions are, however, too long to for the present paper. We specified how h (k) factorizes at 
points -and along lines of high symmetry in the 2D BZ and pointed to the many band crossings and linear 
dispersions ("incipient Dirac cones") caused by the factorizations. Their role, together with that of Fe d- 
As p hybridization for the presence of the d 6 pseudogap at the Fermi level and the details of the shapes and 
masses of the electron and hole pockets were subsequently explained. Thereafter we included interlayer 
coupling, which mainly proceeds via the As z orbitals, and showed how the st and bet 3D band structures 
can be obtained by coupling h (k) and h (k+7rx+7ry) , i.e. by folding the 2D bands into the small BZ. 
This formalism allowed us to explain, for the first time, we believe, the complicated DFT band structures 
of in particular st SmOFeAs, bet BaFe2As2, bet CaFe2As2, and collapsed bet CaFe2As2. What causes 
the complications are the material-dependent level inversions taking place as functions of k z . We found 
several Dirac points near the Fermi level. Whether these points have any physical implications remains to 
be seen. They do not pin the Fermi level, because there are also other FS sheets, and they are not protected, 
because they are merely caused by crystal symmetries, such as the vertical mirrorplane containing the 
nearest-neighbor As atoms, which are easily broken by phonons, impurities a.s.o. 

Copyright line will be provided by the publisher 



42 



O. K. Andersen: On the multi-orbital band structure and itinerant magnetism of iron-based superconductors 



We then studied the generic band structures in the presence of spin spirals, whose Fe moment has a 
constant value, m, but whose orientation spiral along with wavevector q. The formalism simply couples 
ft, (k) to h (k + q) and does not require q to be commensurate. We used the Stoner approximation to SDFT, 
because it is simple and allows one to calculate the spin-spiral band structures as functions of the strength, 
A, of the exchange potential and impose the selfconsisteny condition, A = in (A) I at a later stage. We 
limited ourselves to using this formalism to explain the 2D band structures for stripe order as a function 
of A, but in quite some depth, often using simple analytical theory. What complicates the magnetic band 
structure is the simultaneous presence of As-Fe covalency and Stoner-exchange coupling. That the latter is 
only between like Fe d orbitals, gave the structure of the small-moment SDW-gapping of the paramagnetic 
FS as well as the intermediate-moment propeller-shaped FS . With the goal of eventually understanding 
-and possibly simplifying the calculation of- the spin spiral energy dispersions, we expressed the magnetic 
energy as the difference between magnetic and nonmagnetic band-structure energies, whereby the magnetic 
band structure should be the one corrected for double counting of the exchange interaction. This formalism 
was then applied to the large-moment stripe order and we found its stabilization energy to have two main 
sources (1) the coupling of the paramagnetic xy hole and xy/z electron bands and (2) the coupling of the 
other electron band and the xz part of the doubly degenerate hole band. The Fermi-surface contributions 
to the magnetic energy were found to be comparatively small, and that gave some hope for developing 
a suitable Heisenberg Hamiltonian. We also explained the much discussed coupling in SDFT between 
the stripe moment and the As height above the Fe plane. In the end, we showed and discussed the self- 
consistent spin-spiral moment and energy dispersion obtained from a SDFT calculation, co-authored by 
one of us. EH 
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